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„Je schwieriger das Problem, 
desto größer der Reiz, 
einen Weg zur Lösung zu finden.“ 

Für die allseitig gezogene Ebene mit einem Loche, dessen Rand aus zwei Kreisbögen besteht, wird ein 
neues Lösungsverfahren entwickelt, um das Verhalten der Spannungsfunktion in der Nähe einer Innenecke 
zu untersuchen. Hierzu wird die Scheibe konform auf den Ebenenausschnitt, der durch zwei Gerade begrenzt 
ist, abgebildet. Für die inverse Spannungsfunktion wird ein System von Einzellösungen, die spannungsfreie 
Ränder geben, aufgestellt. Weiter wird gezeigt, wie aus diesen Einzellösungen eine Näherungslösung des 
Problemes zu finden ist, wobei die Randbedingungen eines inneren Schnittrandes und die dort auftretende 
Singularität zu berücksichtigen sind. Ein Glied der Lösung gibt das Verhalten in der Nähe der Innenecke. 
Für das Beispiel mit dem Winkel n/2 in der Innenecke wird das zahlenmäßige Ergebnis angegeben und 

> die Güte der Näherungslösung gezeigt. 

To study the properties of the siress function in the neighbourhood of an inner corner, a new method is 
found for the solution of the problem of the plane provided with a hole, the edge of which consists of twoc acrs 
of eircles, when the plane is under a uniform tensile siress. For that purpose, the dısc is conformally mapped 
onto a part of the plane, bounded by iwo straight lines. For the inverse siress function, a system of particular 
solutions is stated, giving siress-free boundaries. It is shown how to compose an approwimative solution of 
the problem from the particular solutions; for that, the boundary conditions at an inner edge and the occuring 
singularities are to be taken into account. One term of the solution gives Ihe properties close by the inner 
corner. The numerical result is given for the example, that the angle at the inner corner is n|2, and the accu- 
racy of the approzimative solution is illustrated. 

Un plan EtirE uniformement soit pourvu d’un trou, dont le bord est compose de deu® arcs circulaires. 
Pour examiner la fonction de tension pres d’une aröte interieure une methode nouvelle est developpee. A ceite 
fin le disque est repr&sente conform&ment sur un secteur du plan, qui est borne par deux lignes droites. Pour 
la fonetion inverse de tension un sysi&me de solutions particulaires est donne, qui fournirent des bords sans 
‚tension. Il est d&montre, qu’on peut trouver, & l’aide de ceites solutions particulaires, une solution appro=@i- 
mative du probleme en tenant compte des conditions sur une ligne interieure d’intersection et des singularit&s 
y occurrentes. Un des termes de la solution deerit la propriete pres de l’ar&te interieure. La valeur numerique 
est donnee pour l’esample, que l’angle de l’aröte interieure est n/2 la qualite de l’approximation est deEmoniree, 

Paspadarsısaerca HOBBIH MeTon NHecHeNoBaHuA PyHRUHN HANPpAsKeHuA BÖNUSH BHYTPEeHHETO 
yrıa B CIyyae PacTATuBaeMoH BO BCEX HAIPABIeHMAX INIOCKOCTH C OTBEPCTHeM, OTPaHHyeH- 
HBIM ABYMA AyTaMmu epecekamImnxcA OKPY:KHOCTeH. Ü 9TOH MeNIBIO IMIACTHHA KOHhOPMHO 
OTOÖPpasKkaeTcaA HA YACTB INIOCKOCTHU, OTPAHHYEHHOH AByMA NIpAMEIMmm. lust OTOÖpakeHHoü 
PyHXIHH HampsskeHunf MaeTcH CHCTEMA YACTHEIX PemeHHMÄ, YAOBIETBOPAIOIMHX YC/IOBHO 
CBOOONHEIX KpaeB. 3ATeM IIOKA3BIBAeTCH, KAK U3 ITUX YACTHEIX PeIHeHHH MONKHO ÖBITB HaH- 

R) AeHO HPHÖNM?KEHHOE pemieHme MaHHOH 3amayn ıpu yyeTe KpaeBbIX YCAIOBHÄ Y BHYTPeHHero 

yrTHa MH uMemwmeich Tam cuhryaapHoctk. OnuH WıIeH pemeHnA BEIPakaeT TIOBeneHne 
byHEIuu BÖNNSU BHYTpeHHero yrıa. la cıyyası, KorAa BHyTpeHHnnä yros pasuserca m/2, 
UPHBONATCA UHCHIEHHBIE Pe3YJIBTATEI H ITOKA3bIBAeTCH TOYHOCTb METONA pacyeTa. 


A. Problem ! 4? ! 

Auf eine Ebene wirkt im Unendlichen ein allseitiger 

gleichmäßiger Zug p. Die Ebene hat nach Bild 1 ein Loch, 

dessen Rand aus zwei Kreisbögen besteht, so daß Innenecken 

entstehen. Die Spannungsverteilung in der Nähe der Innen- 

ecke soll untersucht werden. Im Beispiel ist der Winkel der 
Innenecke gleich z/2 und die Lochlänge gleich 4 gesetzt?). 


f B. Untersuchung durch Inversion und Verfolgung der 
Er auftretenden Singularitäten 
Zuerst betrachten wir das Verhalten der Airyschen | | | 
Spannungsfunktion #* im Unendlichen. Wir überlegen uns Bildi 


also, welche Glieder F* erhält, wenn wir es für genügend 


1) Es besteht: die Möglichkeit, das Gebiet der Scheibe auf ein Gebiet mit Kreisloch abzubilden und dann 
Lösungen durch Reihenentwicklung zu finden. Dieses Verfahren, das ich auch für die Ebene mit quadratischem 
_ Loch entwickelt habe, wird in einem getrennten Aufsatze erscheinen. Zur Untersuchung des Spannungs- 

zustandes in der Umgebung der Innenecken ist es nicht geeignet. 
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5 pr*?, die weiteren Glieder müssen so beschaffen sein, daß ihre zweiten Ableitungen mit 


wachsender Entfernung gegen Null gehen. Wie beim Kreisloch kommt folglich das Glied 
— cin r* hinzu und weitere Glieder von der Gestalt Re 2*"" und Re r*2*”, wobei die Expo- 
nenten —n negative gerade Werte annehmen. Wir erhalten hiermit. 


Pr= Ss prr—elnrtt A A FEN EAN I re (1). 


Den Wert e wird uns unsere Untersuchung geben. (Ersetzen wir unser Loch durch einen Kreis 
vom Halbmesser 1,2, so wird c=0,72p.) Dieser Ansatz gilt nur bis zum Kreise durch die 
Ecken A und (, also nicht bis zum Lochrande. Innerhalb dieses Kreises müssen wir F* auf 
andere Weise ausdrücken. Am unbelasteten Rande wird dF*—=0, auch F* setzen wir dort 
gleich Null. 

Bequemer werden die Betrachtungen, wenn wir nicht das Gebiet der gelochten Ebene, 
sondern ein geschlossenes Gebiet zu untersuchen haben. Wir bilden darum das Gebiet der 
z*-Ebene konform durch Inversion auf das Gebiet der z*-Ebene, Bild 2, ab, indem. wir 


ul 1 
NY ee nah 
23 = Tier yre A 
setzen. e 
Aus der Spannungsfunktion F* wird für das Gebiet der «y-Ebene eine neue Spannungs- se 
funktion F gebildet. Wir schreiben hierzu F* = r*? o,+9,; hierin sind 9, und @, Potential? , | 
funktionen. Dann setzen wir F —= F*/r*? und erhalten: 


F=o,+r?9. A 
Da bei der Inversion Potentialfunktionen wieder in Potentialfunktionen übergehen, so ist F 
ebenfalls eine Spannungsfunktion. Aus Gl. (1) erhalten wir 


F=Sp+erirt: TERN nn YERSE 


Am Rande wird F und d F gleich Null, so daß auf die Scheibe nach Bild 2 keine äußere Belastung 
wirkt. Dafür tritt imi Punkte U mit &=y=0 eine Singularität auf, deren physikalische Be- 
deutung wir untersuchen. Hierzu berechnen wir die Verschiebung uin z-Richtung. Die Quer- 

zahl r setzen wir gleich Null, da sie auf die Spannungsverteilung unseres 
' Problems keinen Einfluß-hat. Wir erhalten aus Gl. (2): 


2Gu=6 Pen + RIEMEN E 


Für 2=0 also: 
BE EIN N RE Le 


Die zwei Linien UB und UD.des Bildes 2, die bei der unbelasteten i 
Scheibe auf einer Geraden liegen, bilden infolge der Verschiebung « l 
den Winkel 27zc/@. Dent Spannungszustand des Bildes 2 erhalten wir > 
also, indem wir entlang den Linien U Bund UD dünne keilförmige Strei- 
fen herausschneiden und die rechte und linke Hälfte der Scheibe wieder 
verbinden. 

Stellen wir für die rechte Scheibenhälfte eine Lösung auf, und 
nehmen wir für dielinke Hälfte dieselbe Spannungsfunktion, bei der nur 
x durch — x ersetzt ist, so stimmt am Schnittrande die Normalspan- 


BE nung überein. Wir müssen weiter verlangen, daß. am Schnittrande . 
die Schubspannung r7,, = nn... re — 0 oder auch 2 _ — ( ist. Für den Ansatz nach Gl. @) Y 


en das zu; wir bilden für das Glied er?In r Br Ableitung nach R und erhalten are 
Be (er: Inr)=ec(2«21nr-+e). Hieraus folgt, daß für die Linien UB und un, — „ern ") geich RR ” 


Na wird. - 
Wir stellen fest: Entlang der Linien UB und UD wird bei a. Scheibe nach 1 Bi 24 


2 Gun = — 2re|y) RR RER 


ee) 
09%) 2-0 


und 
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Untersuchen wir die rechte Hälfte der Scheibe, so geben uns Gl. (4) und (5) die Randbedin- 
gungen für den Schnittrand DU B. Wichtig ist hierbei, daß beide Bedingungen dimensionsgleich 


sind, denn man findet es durch eine Differentiation von F nach x, 2@u durch zwei Differen- 


tiationen nach % und einer darauffolgenden Integration nach «. 
An Stelle der Gl. (4) ist es zweckmäßig, die Summe der Gl. (4) und (5) zu nehmen: 


oF 
2 @u.0+ (2) SE LER SR RAR EEE EN: 
a) 
Dieses folgt aus der physikalischen Bedeutung der Ausdrücke der linken Seiten. Bilden wir die 


0] { (] 
gleichwertigen Ausdrücke —@ BE FIIRE = und @ n he und integrieren wir beide nach y, 


so folgt: a8 RN 
Br ; } F der a 
-| (+ 2.) dr = —— - [a may=- | day | TE 
0) 0 0 
und 
| a 2 aa oF 
| u: 0% u v 
- ER = Be Ne 9 —E 
lei au [62 „ir et a) EINER 
> Die erhaltenen Integrale sind aber die linken Seiten der Gl. (5) und (6). 


Nunmehr wenden wir uns wieder der Funktion F des Bildes 2 zu. Wir wollen den Span- 
nungszustand im Punkte A der Scheibe des Bildes 1 oder der Scheibe des Bildes 2 untersuchen. 
Wir wissen, daß in Innenecken die Spannungen gegen Unendlich gehen. Wir müssen also die 

. Innenecke stark vergrößert annehmen und die dort möglichen Singularitäten untersuchen. Bei 
dieser Vergrößerung geht das Eckgebiet in eine Scheibe über, die durch zwei Gerade begrenzt ist. 
Um diese Untersuchung genauer durchzuführen, verwandeln wir das Gebiet der Scheibe des 
Bildes 2 durch eine weitere Inversion in den Ausschnitt einer Vollebene, der durch zwei Gerade 


IC 
begrenzt ist. Hierzu gehen wir zuerst zu den Koordinaten ?—=2-F 53 mit dem Anfangspunkt ÜC 


über. Dann führen wir für diesen Punkt die Inversion durch mittels der Umformung 

a — 2 1° Unser Gebiet geht in den Ausschnitt der Vollebene nach Bild 3 über. Wir 
2 er, m 

führen noch e Koordinaten © mit dem Anfangspunkt in der 

Spitze des Ausschnittes ein und erhalten: 


1 
_——2 
das 1-2 (7) 
oder 
Rrer 
SIR EEE ar Ba 5 


In den Punkten U und 0 treten 'Singularitäten infolge des Gliedes r?Inr auf. Wir wollen sie 
hier nicht näher untersuchen, sondern nur ihr Vorhandensein feststellen ?). 
i e Aus F finden wir durch die vorgenommene Inversion für den Ausschnitt der Yollebene die 
neue Spannuhgsfunktion 
Each 
= oder SD N. 6. (9), 


2 


Diese Spannungsfunktion F, und ihre Ableitung oF,/on werden am Rande gleich Null. In 
a inuren o= e* und « erhalten wir die Randbedingungen: 


oF, 

(Fa)uma = 0 und Be Ä) ET 5 10) © 
ou /n=a 

E 2) Bi meinem Aufsatze Rs» mit Kreikbogenkerbeit Z. angew. Math. Mech. 1940 entwickelte ich 

ein Verfahren, bei dem diese Singularitäten und weitere Singularitäten außerhalb des Scheibengebietes unter- 

b werden. Dieses Verfahren könnte auch für das vorliegende Problem angewandt werden. Es würde aber 

Auskunft über das Verhalten der Spannungsfunktion in der Umgebung des Punktes A geben. 
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Wir werden im Abschnitt CO für den Ebenen-Ausschnitt ein vollständiges System von Einzel- 
lösungen angeben, die diese Randbedingungen erfüllen. Diese Einzellösungen multiplizieren wir 
mit noch frei wählbaren Beiwerten. Diese Darstellung der Lösung stellt eine Entwicklung der 
Funktion F, in eine Reihe nach o dar, wobei die Exponenten von o nicht ganzzahlig zu sein 
brauchen. Bei einer beschränkten Anzahl von Gliedern erhalten wir nur eine Näherungslösung. 
Dann gehen wir wieder zur Spannungsfunktion F über und bestimmen die Beiwerte so, daß die 
Randbedingungen Gl. (5) und (6) näherungsweise erfüllt sind. Ein bestimmtes Glied der Summe 
gibt uns dann das Verhalten der Spannungen in der Nähe des Punktes A. Aus dem Wert der 
Spannungsfunktion F im Punkte U erhalten wir weiter den Wert c in Abhängigkeit von p. 


0. Ebenen-Aussehnitt mit spannungsfreien Rändern 
a 


Bild 4a bis d zeigt den Ebenen-Ausschnitt für «—= UNS Zi yR und z. Sind die Rän- 
oF 
der u—= + spannungsfrei, so ist für diese F, = a Nunmehr bilden wir Einzel- 


lösungen, die diese Randbedingungen erfüllen. Hierbei beschränken wir uns auf zur &-Achse sym- 
metrische Lösungen. Wir setzen mit &=&+ in = + und =E—in= e*—-i# für eine 
Einzellösung: 


Er 2 [a ({r+1 + {r+1) + 02 (er—1 + Zr1)] ER Er RT ® 
Beide Glieder in der eckigen Klammer sind in e von gleicher Potenz; a und n sind so zu wählen, € 
daß die zwei Randbedingungen erfüllt sind. In Polarkoordinaten lautet der Ansatz der Einzel- 
ne: F,=0o"t![acos(n+1)a+cosmn—1)ul :.. ...... (12). 
Die Werte «a und n werden aus den Randbedingungen für u—= + berechnet. | 

Aus (F,)u=a —=0 folgt: acos(n+1)a + cos(n — 1) =0, 


aus 5) —=0 folgt: a(n+ 1)asin(a+1)a+(n— 1) asin (n— 1)a=0 (13). 
u=& 
2) Par 
$ 
s b) 
3 
c) 
d) As 
3 € 
DIE Bild 5 


Falls n gefunden ist, so erhält man aus der ersten Gleichung. 
_eos(n—1)a 
 cos(n+1)a 

Um n zu finden, beseitigen wir a aus beiden Gl. (13) und erhalten nach Umformung: 
2na Ka: 2& 
sin2n« sin2«a 


hl, 9, RUE he SZ ee RE ER 


* > ' je 
Beta A ME ie Be PT Zr TA 


ae 
e AA het 
2 a 
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Die rechte Seite gibt die Konstante C = — ‚ die von a, dem gewählten halben Aus- 


sin 2«& 
schnittwinkel abhängt. Diese Abhängigkeit im Bereiche O<a <” zeigt Bild 5. Für «=0 
e . ER j 
wird C= —1, füra= 5 wird C= ©, dann nimmt C bis etwa 4,61 für «—=4,5 ab, um 


wieder bis O= + ® für «= anzusteigen. Für einen festliegenden Wert a, bzw. C sind 
nach Gl. (15) die Lösungen für 2 n«& zu finden. G]. (15) hat außer reellen Lösungen auch kom- 
plexe von der Form 2(n’ + in"). Diese komplexen Lösungen müssen wir beachten, damit 
wir ein vollständiges System von Lösungen erhalten. 


Wir setzen n=n’ + in" in Gl. (15) ein: 
2 (n +in’)a=(sin2(n’ +in')«. 


Die Gleichung zerlegen wir in Real- 
und Imaginärteil; dieses gibt nach 
Umformung: 


F 2n'a 
cos2na= Conawa (16) 
und 
MERAN & 


Durch beide Gleichungen ist 2n”& 
als Funktion von 2n’a gegeben. In 
einem Koordinatensystem der Ver- 
änderlichen 2n’& und 2n"« tragen 
wir beide Kurvenscharen auf, die aus 
den Gl.(16) und (17) folgen. Zur 
ersten Kurvenschar gehört auch die 
horizontale Achse 2n"&=0, zur zwei- 0:4 | 
ten die vertikale Achse 2na—=0. In er‘ 
den Bildern 6a bis d sind die Kurven 


für &1= 0, C, = 1, &o = 1225, 
Gt 59, = 5 oder m, (, =, 
a cams 2,3065, (4,712 "auf- : 


gezeichnet. Die Kurvenschar nach 
Gl. (16) ist ausgezogen, die Kurven- 
schar nach Gl. (17) gestrichelt. Die 
Bilder zeigen deutlicher als Worte, 
wie sich die Kurven mit & bzw. C 
ändern. Die Schnittpunkte beider aulea 

Kurvenscharen geben die gesuchten 

Werte 2na=2na-+:2n'a. Der Wert 2n’«a--2n"a—=0 hat keine Bedeutung, da hieraus 
keine Lösung für F, folgt. 


Um die Vollständigkeit des Systems zu zeigen, ist die Differentialgleichung in den Koordi- 
naten A= In 0, u auszudrücken und durch ein System von Differenzengleichungen anzunähern. 
Für letzteres läßt sich zeigen, daß man mit einem Ansatze, der unserem entspricht, alle vor- 
gegebenen Randwerte und Randneigungen für A—=4, und A= 4, erhalten kann. Für das System 
der Differenzengleichungen ist das System vollständig, also auch für die durch Grenzübergang 
entstehende Differentialgleichung. 


Für komplexe Werte für n werden auch die dazugehörigen Werte für « komplex, aus. 
Gl. (14) folgt: 
cos (n’ + in" —1)a 
cos (n + in” +1)a (18) 
eos (n" —1)a&ojn" a — isin (n’ — 1) Sinn” a Na ; 
eos (m +1) Cojn” a —isin (n + 1)a Sinn’ a 


a=a-+ia = 


198 Weber, Allseitig gezogene bee mit Zweibogenloch Ba. 31 air Er 


Durch Trennung in Real- und Imaginärteil wird daraus «’ und «” gefunden. Für unser Beispiel 


Mit !A.== Im lauten die ersten drei Lösungen, die wir mit den Indizes I, II, III NEXHBDENE 


In a= 0,8177 — 2,562 — 24: 0,544 dr =0,5413, “ 
Inır &—= 2,4460 440,350 = 7,684451,10—20(1,632+4 50,234) ar = — 0,0639 4 50,503 
nr a= 4,4610 + 1 0,5636 = 14,015 451,77 = 2a (2,9784 50,375) ar = — 0,0299 4 50,706. 


BT, Ist n und damit a komplex, so stellt sowohl der Realteil als auch der Imaginärteil des Aus- 
druckes nach Gl. (11) eine Lösung dar, so daß jedem komplexen a n zwei Lösungen für F, 


entsprechen. Zur Berechnung dieser Ausdrücke zerlegen wir — se: +{”-1) in Real- und 
RE 


= Ma + @n-1) — 
1 


I. 3 : ee ee m ee (oe) u a, 


| 
[e=*" # (cos (n — 1) u + isin (m — 1)u) + e*"* (cos (n’ —1)u —%sin (n —1)u] 
— or — [cos nA Cof n” u cos (n’ — 1)u + sinn" A Sin n"usin (n’ — 1) u] 4 
+ vor —1 [— cos n" A Sin n’ usin (n —1)u + sin n"A Eoj n"u cos (n — 1). 


Sowohl Realteil als auch Imaginärteil sind entsprechend unserm Ansatze symmetrisch zur 
&-Achse. 


In gleicher Weise wird — est +£"*") zerlegt. Weiter setzen wir a—= a’ + ia" und er- 


halten den Realteil von F;: 


Fin = Pıt 0°9s- 
Hierin ist: 
91 = a ge" H1 [cos n" A Cof n"u cos (n’ + 1)u + sin n’A Sin n’usin (n +1) 
+ a" g" +1 [cos n"A Sinn" usin (n’ + 1)u — sin n’ A Cof n"u cos (n’ +1) (19) 


und 
9 = 0" 1 [cos n" A Coj n"u cos (n + 1)u + sin n" A Sin n’u sin (n’ — 1) u] 


Wir können g, und 9, auch als Realteile komplexer Funktionen va (A + iu) und fs, (A+ iu) | 
betrachten. Dann lauten die konjugierten Potentialfunktionen y, und %,, das sind die ul 
närteile derselben Funktionen f, und f,, wie folgt: 


y,—= a 0" +! [cos n" A Cof n"usin (n + 1)u — sinn" A Sin n"u cos (n’—+ ul 7 
+ a" oe +! [—cosn"A Sinn" u cos (m + 1)u — sin n" A Coj n"usin (n’ + a, (20) 
pp = or [cos n" A Cof mn” „sin (n’ — 1)u — sin n" A Sinn" weos (m — 1a] - Fi 
Ich gebe die Ausdrücke an, da sie später gebraucht werden. ; . u 
Ebenso schreiben wir für den Imaginärteil von Fu: Bm: ‘ 
_ Fim=Yı+ 0°92. BEE Er 
Hierin ist: A 2 
91 = ad er Hl [— cos n"A Ein n"usin (n" +1)u+- sin nA nenn (1) 
—+ a” 0” +1 [cos n" A Eof n" u cos I +1)a + sinn" Ein n"usin (m Ko 


und 
= or 1 [— cos n"A Sinn” u sin (n’— BL n Ye n" u cos (m’ 1). 
’ - ’ ag n% Ei s 
y,— a og" tr! [cos n" A Sin n" u cos (m +1)w-+sin m" Eoi ei (m + 1) 
+ a" o"*! [cos n" A Coj n'usin (n’ + 1)u —sin n" A Sin n"we cos (m ey ! 

Y, = 0"! [cos m"‘ Sin n"u cos (m — 1)u-+sin n"h Sof "usin ( OR a 


Wir werden die Ausdrücke später für un Wert e=L [ 
sich alle Ausdrücke, da dann cos MA und sin IE 


Die konjugierten Potentialfunktionen sind jetzt: 


Eu a 
E 7 R nn . a (a4) re iu 
y gi r Sour y) sn > JRR 
ee = x, 5 “ E « 2 
} \E . n 
I e 


ı v.. 
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ws Sowohl für die Realteile als auch für die Imaginärteile numerieren wir laufend die Funk- 


eoden 9, und 9, und bilden mit noch freiwählbarem Koeffizienten — ° 5 00% die Potential- 


funktionen 
. IT I 
= 3 2 ‚= a: Pır und d,=—,0 92 % Par: 


_ Hierbei sind in beiden SER die Werte a, dieselben. Aus Di und 2, wird jetzt die Spannungs- 
funktion für F, für den Ebenen-Ausschnitt gebildet: 


3.0 | F,=0,+0:0,. 

Für das Gebiet o<1 sind hierbei für n’--1 positive Werte zu nehmen, für das Gebiet o>1 

.,.. negative Werte. In dem Schnittrande oe=1 müssen beide Funktionen ineinander übergehen 
und die Singularität an der Stelle o=1, «0 bilden. Dieses werden wir jedoch in der & y-Ebene 

durchführen. Hierzu bilden wir aus F, die Spannungsfunktion F der Scheibe in der &y-Ebene. 


Ns; 


7-9) 0, + (2-2 +1) BD: ?=r?d, + (r® 2 +nR.| 
| (23). 
1 19 " 
Pe tetllatlemer landen | 


Hierin sind r’ und r" die Abstände von den Eckpunkten C und 4. 


Dar - Wr ln ae 


Aus F finden wir FEN und 2@u,-0 + e) ‚um diese in die Randbedingungen des 
0 \ 


’ 0-0 
| Schnittrandes Gl. .$) und (6) einzusetzen: 
b oF 1\,00, 08, od oD, 
E20, +0)+—®) la a + ee) a2 29) 


Für &=0 folglich: 


En ee. 
Weiter ist mit 9 = 0: j x 
alle aan arme el 
RE An 
2 N 2 29 (rt + TE in, RE ( 
ee ö Hierin jet a eine. ne: Integrationskonstante, die einer arten Verschiebung ent- 


s Spricht, Für er 0 Ba j 


Y es = 2 (#9, — Do)eo Ay t Folo=n — = ER 
= &=0 


Zeit e # Be; a (4) 


v0 


a 7 c 
h a 5m = —Azelyl DAN a: (25). 
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Aus diesen beiden Gleichungen sind die Beiwerte der Einzellösungen und a, zu finden. Nehmen 
wir eine Näherungslösung mit den zwei Einzellösungen %; und N, njr 4 inzr, so stehen uns 
drei Beiwerte und a, zur Verfügung. Wir können folglich auf der Linie U B die beiden Gleichungen 
an zwei beliebigen Stellen befriedigen. Dieses ist einfacher als die Anwendung eines Minimal- 
problems, da wir damit die analytisch nicht durchführbaren Integrationen vermeiden. Zur 
Erleichterung der Rechnung gehen wir hierzu wieder in das Z-Koordinatensystem über und j 
drücken die beiden Gleichungen in oe und u aus. Die Punkte x=0 entsprechen den Punkten | 
11-—£ 


oe=1, weiter ist nach Gl. (8): z= irre für 0, o=1 folglich | 


Um ni zu finden, überlegen wir uns, daß die Richtung von dx in die Richtung von —do 


übergeht und infolge der Inversion 


do LAD, 1 
ENDET 1 
he 
ist. vr 4 
Die Gl. (24) und (25) lauten hiermit ( 
08, , 908 Bu 
(B, — D,),=1 ug | d6 + =) == 0 NER NE (27) 
und 
2(®,+ P.)=1 5 cu+ 2185 (Fi +P)e1=—5 c2 ig, NE Free (28). 


Auswertung für das Beispiel und Ergebnisse 


1 i 
Gewählt sind für o=1 die Punkte vu = 4° und = a; die Beiwerte a, und a, sind so 


4 

berechnet, daß an diesen - Stellen 

Gl. (27) und (28) erfüllt sind. Für 

andere Punkte des Schnittrandes 
weichen dann die Werte der linken 

Seiten der Gleichungen von den Wer-: 

ten der rechten Seiten ab. Das Er- - 

2 Mi gebnis ist im Bilde 7a und 7b dar- 
zn. (9 ).| gestellt, hierbei zeigt Bild 7a den 


tg 3 und die erhal- l 
5 
| 
4 


Näherung hierfür Verlauf von 2 


)) tene Näherungslösung hierfür. Die 

Differenz stellt den Fehler dar, der 

‘ in einem fünffachen Maßstabe im 

u Bilde 7b aufgetragen ist; ebenso ist 

i A ” % % RE 38 in diesem Bilde die Abweichung der 

linken Seite der Gl. (27) vom Soll- 
werte Null aufgezeichnet. 


mil öf Für die Sehe eu gie E- 
Fehler für elföter (E] . den wir F: ya EN 
x:0 Br 


2 
F.- y=0 — 2 (Bd, + BaE, 
0 b) E 
% 7,4, 1 (,48,.-ı un=—, C 0,425. i% 2 
an ee 
5 5 Da der Wert für F an dieser sa 


gleich — ! „pP I so. Bet 


Bild 7au.b 


2; 
iR 
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Die Spannung bei Punkt A des Ebenen-Ausschnittes wird unendlich und zwar gilt in der Nähe 
dieses Punktes 
u Pa 5 c 0,545 0'544 (— 0,5413 cos 1,544 u + cos (— 0,456) u). 

Da der Exponent von o kleiner als 2 ist, so gehen die zweiten Ableitungen im Punkte o gegen & 
wie 09456, Die Singularität der Funktion F überträgt sich bei der Inversion auf die Scheiben- 
probleme der anderen Gebiete, der Exponent von o und die Abhängigkeit von u bleiben erhalten, 
nur der Zahlenwert der Funktion ändert sich, falls der Abstand vom Umkehrpole nicht gleich 1 ist. 

Für Punkt A des Bildes 2 erhalten wir denselben Ausdruck für F, wenn wir dort o vom 
Punkte A und « von der nach links gerichteten Achse abmessen. 

Für den Punkt A des Bildes 1 wird 


3 1 \1,088 
De = 5° 0,545 nn 0154 (— 0,5413 cos 1,544 u + cos 0,456 u). 
tel 
Hierin ist & der Abstand des Punktes A vom Umkehrpole im Bilde 2 und der Exponent 1,088 


ist der Wert für 2 n,. Hiermit ist das Problem gelöst. 

Neuartig ist bei diesem Verfahren das Aufstellen der Einzellösungen für den Ebenen- 
Ausschnitt und die Berücksichtigung der Randbedingungen eines inneren Schnittrandes und 
der dort auftretenden Singularität. 


Eingegangen: 15. 6. 1950. 


Graphische Integrationsverfahren für die Meißnerschen 
Linienbilder bei Differentialgleichungssystemen und 


Differentialgleichungen n=ter Ordnung 
Von Rudolf Ludwig*) in Braunschweig 


Für Differentialgleichungssysteme und Differentialgleichungen n-ter Ordnung werden graphische Inte- 
grationsverfahren angegeben, die durch Erweiterung der Methoden entstehen, die in früheren Arbeiten ent- 
wickelt worden sind. Die Fehlerabschätzung ist einfach durchzuführen. 


For systems of differential equations and differential equations of the n-th order, graphical methods of 
integration are given, that are generalizations of the methods explained in previous papers. The estimation 
of the error is easily accomplished. 


L’auteur donne des methodes graphiques pour l’integration de sysiömes d’Equations differentielles et 
d’Equations differentielles d’ordre n; elles sont developpe£es par generalisation des meihodes publides precedem- 
ment. L’&valuation d’erreurs se fait simplement. 

 Aasortcs rpaßuyeckne Meroxst uATeTpupoBauna cucrem MahepeunnansHex y paBHeHuH u 
AaudepeHnnanpHbIX YPpaBHeHnf mopanka H. TH MeToAbI ABIMIOTCH PAsBHTHeM METONOB, 
paspaÖoTaHHEIX B IIPeHCHuUX pa6oTax aBTopa. Omenka IMOrPeImHOCTH MOHCET ÖHITB IPOCTO 
mponsBe7eHa. 

Einleitung 
In zwei Arbeiten!) wurden neue Methoden entwickelt, um mit den Meißnerschen Linien- 
bildern bestimmte Integrale und DGln.?) 1. Ordnung graphisch zu integrieren. Die Verfahren 
für DGln. 2. und höherer Ordnung gehen aus denen für DGlIn. 1. Ordnung hervor, während 
Meißner (I [1]) einen anderen Weg wählt, der sich aber schwieriger analytisch erfassen läßt und 
für den auch keine Fehlerabschätzung existiert. Für die hier dargestellten Verfahren ist eine 


‚solche durchgeführt worden in Anlehnung an eine von Tollmien (II [2]) stammende Methode. 


1. Allgemeines zur näherungsweisen Integration und Restabschätzung von 
Differentialgleichungssystemen 1. Ordnung 
Gegeben sei das DGlIn.-System 


d f 
Teeh@n. .. Ya). (k=1, ... n) uf then vige WLzeBine (1) 


*) Aus dem Institut für Angewandte Mathematik und Darstellende Geometrie der Technischen Hoch- 
schule Braunschweig (Prof. Dr. F. Rehbock). 
1) „Analytische Untersuchungen und konstruktive Erweiterungen der graphischen Integrationsverfahren 
von Grammel und Meißner“, Z. angew. Math. Mech. 31 (1951), S. 120—130 und „Graphische Inte- 
ationsverfahren für die Meißnerschen Linienbilder bei Differentialgleichungen 1. Ordnung“, Z. angew. 
th. Mech. 31 (1951), S. 169—178. Hinweis auf diese beiden Arbeiten werden durch Vorsetzen einer „I“ 
oder „LI“ gekennzeichnet, | 
2) DGl(n) = Differentialgleichung(en), Gl{n) = Gleichung(en). 


mit den Anfangsbedingungen 


a CR PR A FAT E 
der n Variablen Yj, - .., Yn. Von den n Funktionen F „ werde RL daß sie in einem Gebiet 
G|a— | <a, |yr— Yu | Sb. 1. beschränkt seien: 


PT.» DS Re Re x,.(2) 
und 2. daß für sie die Lipschitzbedingungen gelten: 


fr (x, Yı--» y)— Fr (8 Yı» ... In) z M 2 Yı — pl; (= 1% she, n) er a Re (3). 


Setzt man die Differenzierbarkeit der Funktionen f,, voraus, so gelten die einfacheren Bedingungen 


of P 
<M, KALT 3), 
132 ) @) 

Die exakten Lösungen des Systems (1) mit den Anfangswerten (1*) lauten: 

8 
Yk (%) I YRs Fe % Yı (8); In (2))d® N I RR (4). 

Dem System (1) soll das Näherungssystem (s. auch 1 [3]) | 
d 
May... Wfl de De c 

en werden, dessen Lösungen den geänderten Anfangsbedingungen | 

Yy;(&) = Yo Yro + Ör a ne Me F ES (6) | 
genügen sollen. Die Lösungen dieses Systems seien y;,(&), es gilt also: 
* ee 
(2) — 90 [I (8 yıl®), -- »» Inla))dr | | 
% 


= /f(® ld), 2. mla))dr er d(&, Filz) +. n(2))dx 


Dieselben Funktionen J;(z) sollen als (nullte) Näherungswerte zur iterativen Lösung des ge- | 
gebenen Systems (1) gewählt werden. Es gilt daher: 


2 _ | | 
IN — yo mn (©, 31(%); ... In(2))dw er HAT, 3 AT FE NR ur (8). . | 
- %o s : vn Fr 
Aus (7) und (8) folgt: 
© s 
IN (x) — (2) = — I — | lo I1r > M)de. | 
s x %o f 
Für die Integrale wird eine Abschätzung ausgeführt: | 
[9031 230 9) del SO nl ur ee 
%o - 


und man erhält: 35 N 
= OBEN AERO ee. 

Setzt man das Iterationsverfahren fort, so findet man unter Sa Er 23 Lipschitz- BR 

bedingungen: \ “ 


. h Se 


© 
| 91 (=) — ed (@)| < M 12 gu) 2) da 
%p 2 Y e N 
und daraus die Abschätzungsformel: Re MARC. N mie 


Ms 


I) ler I re er EHE ER K.. 
Pr 1 
Erfüllen die Näherungslösungen 7; die gleichen Antangsedingungen wie dis: exakte Las, abe er 
also ö.—0, so gilt: ER AB, 2 ee 


Wie in (I) und (II) werde vorausgesetzt, daß die Werte der Funke 
. (geometrisch konstruktiv Auer Be‘ einer SEHR von der Be in 
u 5% oh 


wi 


‚also 


soll das Verfahren gezeigt werden. Die Konstruktionen bieten prinzipiell nichts Neues. Während 


_ man genau so wie bei den einzelnen DGlIn. Integriert man graphisch also z.B. nach dem Sehnen- 


ew. Math, NR 03 
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© gewonnen werden, dabei bedeutet = (mAh), Die Lösung des Näherungssystems kann 


man schreiben: e = 
mtih_ r 

Vk, m+1 — Yk,m = a J(®, Yıla), . es 9()) dx zz 

Im -& 

; &ytih &ytih EN 

= pa. det | Hm I... Nda. ‘4 

I In Er f 

Da alle Funktionswerte nach Einsetzen der Lösung y,(x) nur Funktionen von x sind, kann MB 
man schreiben: BE 
fl I H@)= Fr). Bi 

Aus Se Gl. folgt dann: 1 
Zytih- &yytäh a 


N 9()da—h 3 a, Fulanın) — j F'. (z) )dw. 


Un I 


Dies ist aber wieder der Integrationsfehler RX 
&yrih 


J d,(a)da = = 19% = Ah Or, m» 


1 
Or, m ES, 1% ’ 

wobei RX für die. gefundenen Näherungsfunktionen y,(2) zu bilden ist. Für das gesamte Inte- 

Arsuslersich gilt dann: 

fo, Yu m)da<O! 0<,|RP a, (18 

%o R 


Dieser Wert © ist dann in die Fehlergleichung (11) bzw. (11*) einzusetzen. 


2. Die graphische Integration von Differentialgleichungssystemen 1. Ordnung (Bild 1) 
An dem System von DGln. mit zwei von it abhängigen Variablen « und y 


"=fı(t ©, y)\ 
NEAR u ne a tl) 


bei einer DGI. alle Konstruktionen in einer Ebene durchgeführt wurden, sind hier zwei Ebenen 
notwendig, eine &-.und eine y-Ebene (im allgemeinen Fall » Ebenen y,.). In jeder Ebene verfährt 


verfahren mit Iteration, so bedeutet dies (s. II, Abschn. 3) in ee geschrieben: 
at [fr li io Fr) fı (tert, Ders Ha 5 


Ar— It [fa (te, % 4) + fa (iarı> ur, 9) 5 


Unter Anwendung der Lipschitz- 
en mit der Konstan- 
ten M: | | 

R) of; a 


Ip! 


x-£bene 


sM, h; RR, 


h In = 2) 


Bi; a 
Beralionsyerfahren A: 


Gegeben: «Rh. .heer 
Gesucht: o nach dem Sehnen-Verfahren 


4 
’ 
Ä 
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Bei dem allgemeinen System (1) mit der Lipschitz-Konstanten (3) findet man für die 


Sehnen-Integration: ' 
a) als Konvergenzbedingung 


b) als Fehlerabschätzung: 
h? [23 
= 12 fr (%, Yy +. In) Imaz 
oO 
BR! —__ (en Mle—xı| __ = Eur: 
I Hl< pe BEL 
Je größer also die Zahl der Variablen ist, desto kleiner muß man die Schrittkennzahl A M wählen. 
Bei dem wohl meistens anzuwendenden Sehnen-Tangenten-Verfahren lauten die 
entsprechenden Beziehungen: 


3 3 
a) Migh<,„’ bzw. Mh S, 


(18). © 


h* i 2 Z 2 4 3 A 
9-15 [fe (@, Ye, + 10 fr (%, Yırcı „)-t Ir (2, Yır,+. ER 


b) vr — | Ss Sea 1)5 (Et ee (19). 


Dabei muß der maximale Wert von © unter den » Funktionen f; gesucht werden. 


3. Die graphische Integration von Differentialgleichungen höherer Ordnung 
i Die explizite DGl. n-ter Ordnung 
yn—F(2Y9, 9, Eye 


mit den Anfangsbedingungen 


anach dem Tangenten-Verfahren 


yE (zp) = Yro» 
PR EN 2) 


soll, wie üblich, in ein System von n 
a4. Variablen umgewandelt werden: 


7, I de} y’ Pi ( 
’ 
BL hr 3 
e— ; 
A, (22) 
dyn-1) | 
| IE FERI IE. 
Bild 2, Sehnen-Tangentenverfahren für DGln. 2. Ordnung % 


Die exakten Integrale lauten ganz entsprechend: 


(ya, | Fl yo), yr@))dr 


5 

yo yo | yerdla)dn, (k=0,...,n—2). Be 

ö AR TUR 

Die graphische Integration kann nun genau so erfolgen wie bei den DGl.-Systemen. Da 

aber jede Variable .y%) die Ableitung der vorhergehenden y#-® darstellt, vereinfacht sich die Sc 

Konstruktion, indem man jede der y®-Ebenen mit der vorhergehenden y*-V-Ebene ch 

Drehung um 90° zusammenlegen kann, wie aus den Definitionen der Polaren Gl. I (9) unmittelb 

abzulesen ist. Die Fehlerabschätzungen können unverändert von denen über Systeme von DGln. 
übernommen werden. i Re Ye 
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Für die DGl. 2. Ordnung (Bild 2) 
= f(®, y. y) 


mit den Anfangsbedingungen y(2,)—=y, und y’(2,)=y, soll das vorteilhafteste Sehnen- 
Tangenten-Verfahren noch einmal für die einzelnen Schritte skizziert werden. 
Das äquivalente DGl.-System lautet: 
a, = 
Fre = Seyy). 


Die Integration soll sich über das Bereich a <x <b erstrecken. Aus den Größen 1, = und u 

Yy 
innerhalb des betrachteten Bereiches hat man die Lipschitz-Konstante zu N: Unter 
Beachtung der Konvergenzbedingung für die Iteration 


Mh< 3/2 
wählt man die Schrittweite A. 


Wie bei den DGlIn. 1. Ordnung (Zusammenstellung II, Abschn. 6) erfolgt die Konstruktion 
in den folgenden Schritten: 


1. Schritt: Bei halber Schrittweite: Konstruktion von Pih und Pıa auf den 
Tangenten der Polaren in P, und Ps. 
Berechnung von yın aus der DGI. 

2. Schritt: PoPif% und P’% Pi) nach dem Sehnenverfahren integriert: Pıy und Pi. 
Berechnung von yj/$ aus der DGI. 
Wiederholung EEE 2. Schrittes bis die Iterationen innerhalb der Zeichen- 
genauigkeit konvergieren. 

3. Schritt: Bei ganzer Schrittweite: Konstruktion von PP und p nach dem 
Tangentenverfahren. 
Berechnung von yı ® aus der DGI. 

4. Schritt: P, Pıya pP, und N 4 Pij pP” nach dem Sehnenverfahren integriert. 
Durch Mittelwertbildung Pi” und P}” bestimmt. 
Berechnung von yı > aus der DGI. 
Wiederholung des 4. Schrittes (Iterationsverfahren). 

5. Schritt: Wie 3. Schritt für das folgende Intervall x, bis x, ergibt PP und PP. 


6. Schritt: Wie 4. Schritt: PY, BERN P, und pP; R p,® BER POFUSW. 


4. Beispiele. 
4,1. Die Differentialgleichung 2. Ordnung für die Kettenlinie (Bild 3) 


Das physikalische Problem, die Gestalt eines der Schwerkraft unterworfenen Seiles zu 
finden, führt auf die DGI. 


y—kyiFy® 
Als Beispiel soll das folgende Anfangswertproblem durchgeführt werden: 
ee Yı+y? ’ 
y0)=1,y()=—1, 
für den Bereich Veran, 


Die Originalzeichnung wurde mit = 0,1, e= 150 mm durchgeführt und ergab folgende 
Werte: 


mer merre I v Ir 
Dr LAD 1,000 —1,000 1,4142 s | 0,8 0,589 0,089 | 1,0040 

0,05 0,936 1,3697 9 | 09 0,581 —0,011 1,0001 
170,1 0,908 0,864 1,3216 10 | 10 0,586 0,113 1,0064 
2 | 02 0,833 —0,739 1,2434 11.1 0,605 0,216 1,0231 
8 1,08 0,759 0,626 1,1798 2 | 12 0,632 0,321 1,0502 
4 |'04 0,707 0,511 1,1230 1821°,8 0,669 0,424 1,0862 
5 | 05 0,656 0,399 1,0767 14 | 14 0,722 0,532 1,1327 
6 | 0,6 0,625 0,293 1,0420 15 | 15 0,780 0,643 1,1889 
7.107 0,599 —0,193 1,0184 


N a a u 
r K rd et mn vier en ne © 


IR ee) 


N 
ur FWITN U Ta Fa) We 
Sm H PN 


=» 

2 
5 
Pe 
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Die Berechnung der Funktionswerte wurde mit einer Quadrattafel durchgeführt. Außer FE 
am Anfang braucht man innerhalb der Zeichengenauigkeit kaum mehr als eine Iteration durch- 
zuführen. Sehr genau läßt sich aus der Zeichnung das Minimum der Funktion ablesen, da dort 
die zweite Polare C, und das Linienbild Z sich schneiden und außerdem Z, dort eine Spitze 
haben muß (I, Satz 5—7). 

Zur Fehlerabschätzung: Wegen f} = y’ und &— Y1-+ y’? ergibt sich für das angegebene 
Bereich die Lipschitzkonstante M=1. Die Konvergenz der Iteration ist gesichert wegen 
Mh <0,1< 3/2. Aus der Differenzentabelle der f, und fs werden die maximalen Differential- 
- quotienten näherungsweise bestimmt: 


Bi \ : KB... == 850, Veen: a 11% ee —r 1. 


03, 0883 


Beispiel 1: 
Kr 
J0)=1,y10)=-7 


Bild 3 


"Mit jet Integrationsfehler | 
h* 173 
I — 10h -H 10 f, -H Ye — 5,5 - 10-% 


schätzt man den Fehler des letzten Funktionswertes ab: 


915 —Fu|<S ern) A hei 


915 = 0,780 + 0,005. 


Der hier leicht angebbare exakte Wert des Integrales lautet: 


y—Y26oj 2 — Sin 2-1 — 1234 also FRI 
der Näherungswert mit den Fehlerschranken schließt also den exakten 
Führt man zum Vergleich dasselbe Beispiel durch unter 
Regel in der Form, wie sie in II, Abschn. 8 angegeben wurde, s 
ausgeführt werden soll: | 22 ER 


also: 


o 
I 0,7834 + 0,0001. 


er {ü “L 
R * 


. wi nr T r b 
ü up ar ee ö a D 
2 ch DR A a ACT IE EET 
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4,2. Die ungedämpfte, anharmonische Pendelschwingung (Bild 4) 


Während für kleine Amplituden die harmonische Schwingung der DGl.: 

e | » (En N en (a) 

mit den Anfangsbedingungen: 
| P0)=p, PO) —0 

‚genügt, wird für größere Amplituden die Schwingung anharmonisch und genügt der DGI.: 

ET 0: De TEE LEN (b). 


K 


[ 


Beispiel 2: 
7244-0490 
J0)7 4 10)-0 


Bild 4 u; 


Er Uns interessiert in diesem Falle die Schwingungsdauer und man erhält aus den Lösungen der 
Ber, . DGln. für eine Viertel-Schwingung: E 
EEE ER a: 

3 ER Fall @): E R i een | und 


ya\ 


’ 
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mit den Anfangsbedingungen: 7 5 
9% — 4 ’ Az 0. 


Ohne auf die Einzelheiten der Konstruktion nochmals einzugehen, erhält man für den 
gesuchten Zeitpunkt des ersten Durchganges durch die Nullage, also für die Zeit einer Viertel- 
schwingung 

7 _ 1,059 
Re 

Diesem Wert stehen gegenüber die Werte aus den DG]. (a) mit T/4—=1,014 und (b) mit 
T]4—=1,055 [3]. Die Vergrößerung der Schwingungsdauer ist also deutlich bemerkbar. Auch 
hier wird die Ablesung der Nullstelle dadurch erleichtert, daß das Linienbild eine Spitze hat 
(I, Satz 5—7). 


Zusammenfassung 


Die letzten Beispiele haben gezeigt, daß.die dargestellten graphischen Methoden geeignet 
sind, um mit ihnen ohne großen Zeichen- und Rechenaufwand schnell eine Lösung einer DGl. 
beliebiger Ordnung mit guter Näherung zu konstruieren. In der dargestellten Weise ist es auch 
möglich unter Verwendung der Meißnerschen Linienbilder sich einen Überblick über die 
Größenordnung des gemachten Fehlers bei sorgfältiger Zeichnung zu verschaffen. Ein daraus 
entwickeltes numerisches Verfahren liefert brauchbare Ergebnisse bei geringem Umfang der 
Rechnungen. & 


Schrifttum 
[1] Fr. A. Willers: Methoden der praktischen Analysis. $18,7. 8.166. Berlin-Leipzig 1928. 
Sind die n Gln. 2, = 9, (®,...%,„) durch Iteration zu lösen, so ist die Konvergenzbedingung: 
nM < 1, wenn IM |<M. 
[2] A. Föppl: Vorlesungen über technische Mechanik. Bd.IV. 10. Aufl. $13, S.86—89. München 1944. 
[3] Jahnke-Emde: Tafeln höherer Funktionen. 4. Aufl. S.53 und 62—67. Leipzig 1948. 


Eingegangen 19. 6. 50. 
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Vollständige Darstellung der turbulenten 


Geschwindigkeitsverteilung in glatten Leitungen 
Von H. Reichardt in Göttingen 
(Aus dem Max Planck-Institut für Strömungsforschung) 


Auf Grund experimenteller Untersuchungen des turbulenten Austausches wurde eine Formel für die 
Geschwindigkeit turbulenter Strömungen in Rohren und Kanälen entwickelt. Diese Formel stellt die Ge- 
schwindigkeit über einem Halbmesser vollständig dar, d. h. sie gilt näherungsweise sowohl für die unmittel- 
bare Wandnähe als auch für den mittleren Teil des Strömungsquerschnittes. 


From experimental investigations of turbulent exchange of momentum a formula is derived for the velo- 
„eity of turbulent flow in tubes and pipes. This formula gives a complete representation of the velocity distri- 
bution, i. e. it holds good close to ihe wall as well as in the core of the flow. 


En vertu d’une examination ewperimentale d’&change turbulant d’impulsion une formule fut d&velopp£e 
pour la vitesse d’Ecoulements turbulents dans de tuyauz et de canals. Cette formule donne une reprösentation 
complete de la distribution de la vitesse, c’est-a-dire, elle est valable non seulement dans la voisinage imme£- 
diate de la paroi, mais encore dans la partie centrale de l’&coulement. 

Ha 0cHOBaHHuH 9KCHEePHMEHTANBHOTO HCCHIeNOBAHHA BUXPEBOTO OOMEHA BEIBONHTCH POPMyıLa 
AA CKOPOCTH BUXPEBOTO IOTOKA B TPyOax u KaHalax. ITa POpMyıa NaeT IONHOe pacıpe- 
AereHne CKOPOCTH, T. €. OHA NeÄCTBUTENBHA HU BONHSH CTEHOR U B CpenHel YACTH IOTORA. 


1. Problemstellung 
. Die bisherigen Untersuchungen der turbulenten Strömungen in Rohren, Kanälen und in 
der Reibungsschicht an der ebenen Platte haben ergeben, daß die Geschwindigkeit annähernd mit ? 
dem Logarithmus des Wandabstandes ansteigt. Die Auswertung dieser Versuche mit Hilfe all- f 
gemeiner Dimensionsbetrachtungen lieferte das turbulente Wandgesetz, das folgendermaßen ; 


lautet: = 
U 


1 ER E. 
u ES en .£ 


Hier bedeuten u die mittlere Strömungsgeschwindigkeit im Abstande y von der Wand, u*—= Yrjo 3 | 
die „‚Schubspannungsgeschwindigkeit‘“, r, die Wandschubspannung, e die Dichte des Mediums, > 


ge 


dı 
0 r 
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n= yu*|v der in Form einer Reynoldsschen Zahl geschriebene dimensionslose Wandabstand 
und »=u/o die kinematische Zähigkeit. x und c sind Konstanten. 

Die Ableitung des Wandgesetzes hat zur Voraussetzung, daß die turbulente Scheinrei- 
bung der Schubspannung 7, in der Nähe der Wand überall praktisch gleichgroß ist. Gl. (1) gilt 
also eigentlich nur für eine dünne Schicht in großer Wandnähe. Diese Schicht darf aber anderer- 
seits auch nicht zu dünn sein, da die Turbulenz in unmittelbarer Wandnähe schon aus Kontinui- 
tätsgründen erlischt und die Strömung dort laminar dahingleitet. Die bisherige Theorie ist also 
von beiden Seiten auf einen schmalen Existenzbereich beschränkt. 

Bekanntlich gilt das Wandgesetz auch noch in größerer Entfernung von der Wand, wo die 
turbulente Scheinreibung 7; von der Wand-Schubspannung r, wesentlich abweicht. Dies ist 
ohne Zweifel von praktischem. Vorteil. Vom theoretischen Standpunkt aus betrachtet ist jedoch 
die Gültigkeit einer Formel außerhalb des Gültigkeitsbereiches ihrer grundlegenden Vorausset- 
zungen solange unbefriedigend, als für diese Tatsache noch keine einwandfreie physikalische Er- 
klärung gegeben worden ist. 

Erst im Innern des Strömungsquerschnittes treten merkliche Abweichungen vom Wand- 
gesetz auf und zwar sind die tatsächlichen Geschwindigkeiten größer als die nach Gl. (1) berech- 
neten Geschwindigkeiten. Wenn diese ‚„Zusatzgeschwindigkeit im Innern‘ (die wir mit u; be- 
zeichnen wollen) auch gering sind, so ist doch die Tatsache von grundsätzlicher Bedeutung, daß 
hier die dimensionslose Geschwindigkeit nicht allein von n, sondern auch vom Wandabstand un- 
mittelbar abhängt. Eine einwandfreie Darstellung dieser besonderen Verhältnisse im Innern der 
Reibungsschicht steht bislang noch aus. 

Der Verfasser hat sich daher die Aufgabe gestellt, die aufgezeigten Lücken der Theorie aus- 
zufüllen. Zunächst soll das Wandgesetz durch die Darstellung der Geschwindigkeitsverteilung 
bis in die Jaminare Wandschicht hinein vervollständigt werden. Dann wollen wir die im Innern 
der Reibungsschicht bestehende Abhängigkeit der Strömungsgeschwindigkeit vom Wand- bzw. 
Mittenabstand klären und formelmäßig beschreiben. 

Wir bschränken unsere Betrachtungen auf den Typ der ausgebildeten Strömung mit Druck- 
abfall, wie er in Rohren oder Kanälen beobachtet wird. 


2. Grundsätzliche Betrachtungen 


Wenn hier von einer vollständigen Darstellung des turbulenten Geschwindigkeitsprofiles 
die Rede ist, so ist damit nicht eine vollständige Lösung des Problems der turbulenten Geschwin- 
digkeitsverteilung gemeint. Exakte Lösungen kann es in der Turbulenzforschung erst dann geben, 
wenn das Turbulenzproblem als solches gelöst ist, wenn also (wie beispielsweise bei den zähen 
Strömungen) die physikalischen Probleme auf rein mathematische Aufgaben zurückführbar sind. 
Von dieser höchsten Entwicklungsstufe ist aber die Turbulenzforschung noch sehr weit entfernt. 


Man kann sich höchstens mit Arbeitshypothesen behelfen, über deren Brauchbarkeit das Ex- 


periment entscheiden muß. i 

Von den bisherigen Bemühungen zur theoretischen Unterbauung der Turbulenzforschung 
sind vor allem die Mischungsweghypothese Prandtls und die Ähnlichkeitshypothese v. Kar- 
mans zu nennen. Was die Rohrströmungen anbetrifft, so führten beide Hypothesen genähert 
zur logarithmischen Geschwindigkeitsverteilung. Wir müssen uns aber darüber im klaren sein, 
daß dieses Wandgesetz keineswegs eine zwingende Konsequenz der diffizilen Überlegungen über 
das Verhalten der Wirbelballen oder über die Ähnlichkeit der Wirbelbewegungen darstellt. Für 
die rechnerischen Ergebnisse war allein die Annahme wesentlich, daß der durch den Ausdruck 

7, [du 

.o/dy 
schien ‚plausibel‘, sie war jedoch keine logische Folge jener Hypothese. Der wissenschaftliche 
Fortschritt lag also nicht in den Hypothesen selbst, sondern in einer plausiblen Annahme über 
einen Hilfsbegriff, die durch Versuche bestätigt werden konnte. 

Für die Weiterbildung der Theorie erschien es uns daher zweckmäßig, irgendwelche speku- 
lativen Betrachtungen zu vermeiden und die für die Ableitungen benötigte grundlegende Vor- 
aussetzung unmittelbar auf die Erfahrung zu stützen. Durch die formelmäßige Wiedergabe der 
zugrundegelegten Meßergebnisse, die ja nicht ohne Phantasie zu bewerkstelligen ist, kommt na- 
turgemäß eine gewisse Ungenauigkeit in die rechnerische Darstellung. Auf jeden Fall ist aber 
ein solcher ungenau wiedergegebener Tatbestand ein besseres Fundament für die Theorie als 
eine rein hypothetische Voraussetzung. 

Für die vorliegenden Untersuchungen wollen wir als wesentlichen Hilfsbegriff die Aus- 
tauschgröße A verwenden, die durch die Gleichung 


du 
} se Sem 


definierte Mischungsweg 1 mit dem Wandabstand linear ansteigt. Diese Annahme er- 


. 2) 


14 


w 


weiteren Rechnungen eingehen. 
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definiert ist (7, = turbulente Schubspannung)!). Die turbulente Austauschgröße ist das formale | 
Analogon zum molekularen Zähigkeitskoeffizienten, dessen Definitionsgleichung 7, = u = ! 
lautet. Die Größe A hat aber nicht nur eine formale, sondern auch eine wesentlich physikalische 
Bedeutung. Denn die Frage nach dem Verhältnis der turbulenten Schubspannung 7; zur zähen 
Schubspannung r,, (also nach dem Verhältnis A/u) hat einen unmittelbar physikalischen Sinn. 
Sie spielt daher bei fast allen Turbulenzuntersuchungen, insbesondere auch bei den vorliegenden 
Betrachtungen, eine wesentliche Rolle. 

Wir beginnen nun mit unseren Untersuchungen und stellen an die Spitze unserer Rech- 
nung folgende identische Beziehung: 


u du 
a) TR 


d N = 2 — ep . . . . - . . . . . * ” . 0 (3) 
dy/o 

(der Index O0 kennzeichnet die an der Wand auftretenden Werte des Geschwindigkeitsgradienten 

bzw. der zähen Schubspannung r,,)?). Wie man sich durch Verifizieren von (3) leicht überzeugen 


kann, bedeutet der Differentialquotient der dimensionslosen Geschwindigkeit nach dem di- 
mensionslosen Wandabstand nichts anderes als das Verhältnis des Geschwindigkeitsgradienten 


‚an der Stelle y zu dem entsprechenden Gradienten an der Wand. 


Wir führen die Gesamtschubspannung r ein: [9 
r -„.t%= Br 3 a Ra RT 


Unter Verwendung dieser Gleichung folgt aus (3) 


a(&) 2 
ER 2° 
Er : He ee ae ee 


Aus (5) ist zu ersehen, daß die dimensionslose Geschwindigkeitsverteilung von den Größen 
° und 
To u 
genden betrachten werden, ist 2 bekannt und zwar gleich 1 — — “ (wenn y der Wandabstand und 


= = abhängt. Bei bo ausgebildeten Strömungen er Druckabfall, die wir im fol- 


r der Halbmesser des Rohres oder Kanales bedeutet). Unsere Aufgabe reduziert sich also auf 
die Frage nach der Abhängigkeit der Austauschgröße vom Wandabstand und von der Rey- 
noldsschen Zahl. 

Wie wir bereits sagten, läßt sich die Austauschgröße ebensowenig aus der praktisch noch 


. unerforschten Turbulenzstruktur berechnen wie etwa der Mischungsweg oder irgendein anderer 


Hilfsbegriff. Wir haben daher die dimensionslose Austauschgröße. A/un, durch Versuche be- 


* a 
stimmt (m = a . Diese experimentell ermittelte Funktion (siehe Bild 5) wird in unsere 


3. Das Wandgesetz 


In Wandnähe läßt sich die Austauschgröße verhältnismäßig leicht aus der gemessenen Ge- 
schwindigkeitsverteilung bestimmen. Bereits Nikuradse?®) hat eine solche Bestimmung von 
A durchgeführt. Seine Messungen ergaben, daß in Wandnähe | 

UN. r 


ist (der Proportionalitätsfaktor zist etwa 0,4). Ingrößerer Entfernung (etwa oberhalb von I= 0 ) 


- 1) Die Anwendung des Mischungswegbegriffes wäre unzweckmäßig, da bei der Darstellung des Vverhät- f 
nisses 7/7 (Ss. u.) der Geschwindigkeitsgradient nicht herausfallen würde. Ferner tritt auch beidenhier 
behandelten Strömungen die bekannte Schwierigkeit auf, daß der Mischungsweg in der Symmetrieache x 
(bzw. -ebene) zu einer Singularität führt. Im mittleren Bereich ist nämlich i in Anbetracht derlinearen Veräufe 


der Schubspannung sowie des Geschwindigkeitsgradienten (s. Bild 4) I= y= du PEN, wenn z den vg 
Abstand von der Mitte bedeutet. 1 wird daher bei z— 0 unendlich. el day” % Vz E 
®) Wir werden fortan den Mittelwertstrich über dem u fortlassen, wenn eine ‚Verwechslung von u Bo Br 
dem Augenbliokswert “= % + w nicht zu befürchten ist. Eh IE 
») J. Nikuradse, Gesetzmäßigkeiten der turbulenten Strömung in glatten Rohren. VITO. ra De 
schungsheft 356, 1932. £ A x 2% I 
a re: ‘E 0 4 S + Par: 
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wächst A/un, schwächer als linear mit y/r. Auch über diesen Verlauf geben die Untersuchungen 


Nikuradses Auskunft (allerdings sind seine Werte oberhalb von I. 0,8 unsicher. [Näheres 
siehe im 4. Kapitel)). 5 

Über den Gültigkeitsbereich von (6) in Richtung zu kleinen Wandabständen ist nichts 
Näheres bekannt. Fest steht nur, daß GI. (6) in allernächster Wandnähe nicht mehr gelten kann, 
weil der turbulente Austausch dort aus Kontinuitätsgründen verschwindet und sich die Flüssig- 
keit parallel zur Wand laminar bewegt. 

Da für das Erlöschen der Turbulenz an der Wand die Kontinuität verantwortlich ist, so 
läßt sich aus der Kontinuitätsgleichung eine Aussage über das Abklingen von A gewinnen. Die 
Kontinuitätsbedingung für die Geschwindigkeitsschwankungen w’ (parallel zur Wand) und v’ 
(senkrecht zur Wand) lautet: 


av ou’ 
| oYy ern m Fr * . . . . . . . . ” . . . . * (7). 
Bezeichnen wir das Verhältnis «’/“ mit &, so können wir für w’ ansetzen: 
; & n 
"= sü=a(7,), BR RR RR 3} 
wobei & eine Funktion der Zeit und der Ortskoordinate & darstellt. Setzen wir (8) in (7) ein, so 
) folgt: 
: a y* 
V a EN HR dy 3 . . . . . . . ... . . * . . (9). 


Während wir also in erster Näherung einen linearen Anstieg von w’ mit % (entsprechend 
dem linearen Anstieg von u mit y) anzunehmen haben, folgt aus der Kontinuitätsbedingung, daß 
v —y® ist. Daraus folgt weiter für die turbulente Scheinreibung 


— 002 (du) 
eure, N EN r 
T; euv' —(, se) s (10) 
Für den Verlauf der Austauschgröße in unmittelbarer Wandnähe erhalten wir dann: 
A 2 *\3 
a (I) DT I RL (ul 
u v 


(7) 


Die Austauschgröße ist also proportional n7? (über die Größe des Proportionalitätsfaktors können 
wir hier allerdings nichts aussagen). Das bedeutet ein außerordentlich schwaches Ansteigen von 
A bei den kleinen n-Werten. Sowohl die erste als auch die zweite Ableitung von A nach n ver- 
schwindet an der Wand. Erst die dritte Ableitung hat einen endlichen Wert *). 

Für die formelmäßige Darstellung von A/u als Funktion von % bzw. n benötigen wir daher 
eine Funktion, deren erste und zweite Ableitung an der Wand Null ist. Andererseits muß diese 
Funktion in einiger Entfernung von der Wand in den linearen Verlauf xn asymptotisch überleiten. 


Ü *) Anmerkung bei der Korrektur. 
An der vorstehenden Ableitung ist zu bemängeln, daß die dritte Schwankungskomponente w’ vernach- 
2 ’2 . 
lässigt wurde. Bei dieser Vereinfachung ist, wie Gl. (11) zeigt, = » n?, wenn = m a #0ist. Diese Bedin- 
gung ist aber nur in nichtausgebildeten Strömungen erfüllt. 4 
Wird die w'-Komponente berücksichtigt und ausgebildete Strömung angenommen, so erbält man 


11): ne! ı De 
statt ( ) . a 1 du’ 10 a 5 ; (11a) 
nn awilan)ae ln), % En 
TR r 
(© = zu */v). Hier ist zu fragen, ob 3) und * ) korreliert sind oder nicht. Falls diese beiden 
on)o 05 \0n)o 


Größen völlig unkorreliert sind, so würde daraus folgen, daß die Austauschgröße nicht mit 7° ansteigt, sondern 
mit einer noch höheren Potenz von n. 
Wir kommen also zu folgendem Ergebnis: Wird der Verlauf der Austauschgröße in unmittelbarer Wand- 
© nähe durch eine Potenzreihe von n dargestellt, so entfallen die Glieder mit 7 und n?. Das Glied dritten Grades 
3 - existiert mit Sicherheit nur bei nichtausgebildeten Strömungen. 

Die Kenntnis des Verlaufes der Austauschgröße in nächster Nähe der Wand ist für die Theorie der Wärme- 
übertragung von wesentlicher Bedeutung. In einer neueren Arbeit des Verfassers über den „Einfluß der wand- 
nahen ande auf den turbulenten Wärmeübergang‘ wurde die Proportionalität der Austauschgröße mit 7? 

_ angenommen (Mitteilungen aus dem Max-Planck-Institut für Strömungsforschung Nr. 3, 1950). Die dort für 
sehr kleine n-Abstände angegebenen kleinen A/4-Werte sind möglicherweise zu groß, da A/u möglicherweise 
mit einer höheren als der 3; Potenz von n ansteigt. Noch wesentlich höhere Beträge der Turbulenzreibung in 
_ Wandnähe setzt K. Elserin seiner Arbeit „Reibungstemperaturfelder in turbulenten Grenzschichten‘‘ voraus 
(Mitt. Inst. Thermodynamik ETH Zürich, H. 8 (1949)). 
AED ? 14* 


u 
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Diesen beiden Anforderungen wird nun die Funktion 


= #ln—mTang 2) a a N 
gerecht. Wir ‘werden daher diese Gleichung als Ansatz für A/u benutzen. 

Die in (12) auftretende Konstante n, ist ein Maß für die ‚Stärke der laminaren Wand- 
schicht“. nı ist so vorzugeben, daß die mit Hilfe von (12) berechneten Geschwindigkeiten der 
turbulenten Wandströmung mit den gemessenen Geschwindigkeiten übereinstimmen. 

Der durch Formel (12) beschriebene Verlauf von A/u ist in Bild 1 dargestellt und zwar 
für x = 0,4 und 7, = 11%). Wie wir noch weiter unten darlegen werden, führt dieser Wert von 
7, zu dem Werte 5,5 der Konstanten c des logarithmischen Wandgesetzes (1). 


75 


SS 


70 
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Bei Verwendung von (12) und unter Beachtung der Tatsache, daß in ausgebildeter Strö- 


mung 7=7T, (1- # ist, folgt aus (5) 


a(&) | 
ni BE Ai U . (3). 
| dn 1+=n— m Tang n/m) 
Die Quadratur dieser Gleichung ist leider nicht mit Hilfe der bekannten Rechenoperationen 
auszuführen. Wir müssen daher nach Aushilfen suchen. 

” “ Zunächst können wir im Zähler y/r gegenüber 1 vernachlässigen, da wir es ja in diesem 
Kapitel ausschließlich mit wandnahen Vorgängen zu tun haben. Ferner empfiehlt es sich, auf 
der rechten Seite von (13) ein Glied 1/(1-+xn) abzuspalten, weil dadurch wenigstens eine Teil- 
integration möglich wird. Nach Ausführung dieser Integration erhalten wir dann 


er 1 j 
.- iz butem)+im RT N 3 (1). 
Hier ist m) durch die Gleichung definiert: 2% 
“mTangnmdn a AR 
RR | (Fan —um Tang la) Ü-Fan) a 


Auch die Quadratur von (15) ist nicht ausführbar. Wir Kuna aber den Integranden f Ve 
näherungsweise durch die Funktion 


"(em 4 (bn—1) ®) wi A 


> ‘7 


4) Wie Bild 1 zeigt, ist A < uim RN n< =11lbzw. A> im Gebiet n 38 # 
lose Wandabstand n, = 11 ist daher ein Maß für die Stärke risk nun 
kulare Reibung die turbulente Reibung überwiegt. 


# 


15 A ei 
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ersetzen, wodurch die Integration möglich wird. Wir erhalten dann für f}: 


H=c (1-emm Bee! en) BEE I Wo EEE Sue 1 FE EZ 
Nı 
Hier sind die Konstanten c, und b aus den Grenzbedingungen zu ermitteln. 
Ersetzen wir in Gl. (14) f durch die in (17) dargestellte Näherungslösung f}, so erhalten 
wir für hohe n (bei denen 1 gegenüber xn zu vernachlässigen ist) 


1 1 
En hnt—,inx#+o A AB TERRA ARENA RN) 7 


ur % 


Vergleichen wir diese Formel mit der für hohe n geltenden Gl. (1) so folgt, daß 
1 
" In#x<+c,=c 


gesetzt werden muß. Unter Verwendung der üblichen Zahlenwerte # = 0,4 und ce=5,5 er- 
halten wir dann 0, =5,5 + 2,3 = 7,8. 

Der Wert der Konstanten n, ergibt sich aus der Bedingung, daß das Integral der Gl. (15) 
bei Verwendung von x = 0,4 für hohe n den Wert 7,8 erreicht. Wir haben Gl. (15) für verschie- 
dene n,-Werte graphisch integriert und festgestellt, daß der geforderte Grenzwert 7,8 erhalten 
wird, wenn man n, = 11 setzt. 

Die Konstante b war aus den Grenzbedingungen der Wandschubspannung zu ermitteln. 
Diese werden nicht wesentlich verletzt, wenn man 5b = 0,33 setzt. 

Führen wir nun die vorstehenden Ergebnisse in Gl. (14) ein, so erhalten wir als Näherungs- 
formel für das Wandgesetz: 


U = er 3 
“= 25m(1404n)+78(1—e are 0,07) RE SA), 


Der Verlauf von w/u* ist für kleinere n-Werte in Bild 1 und für etwas größere n in Abb. 3 
dargestellt (und zwar hier als Funktion von log n). Die ausgezogenen Kurven zeigen das Ergeb- 
nis der graphischen Integration, 
während die gestrichelten Kurven- 
züge den w/u*-Verlauf nach der 
Näherungsformel (18) wiedergeben. 
Wie man sieht, bestehen zwischen 
den durch graphische Integration 
und den durch Näherungsrechnung 
ermittelten dimensionslosen Ge- 
schwindigkeiten an keiner Stelle 
wesentliche Unterschiede. Etwas 
größere Differenzen bestehen stel- 
lenweise zwischen f’ (ausgezogen) 
und f, (gestrichelt), wie man aus 
Bild 2 erkennt. Diese Abweichung 
ist aber insofern bedeutungslos, als Bild 2 
man ja für irgendwelche Rechnun- 
gen, in die die Ableitung von f eingeht, nicht auf die Näherungsformel (16) angewiesen ist, 
sondern den Integranden von Gl. (15) unmittelbar verwenden kann. 

Unsere Gl. (18) hat gegenüber der bisherigen Formel (1) den Vorteil, daß sie die turbu- 
lente Scheinreibung für jeden n-Wert des wandnahen Bereiches beschreibt. Während das u/u* 
der Gl. (1) mit abnehmendem n nach —co geht, ist unsere Formel auch bei n = 0 richtig, 
indem sie hier die Haftbedingung u = 0 erfüllt. 

Für sehr kleine endliche Werte von n folgt aus (18) 


w—n 


"Dies ist der lineare Geschwindigkeitsanstieg in der laminaren Wandschicht, wie wir ihn auch 
unmittelbar aus Gl. (5) erhalten, wenn wir dort A= 0 setzen. Dieser lineare Anstieg reicht 
bisn #4 (s. Bild 1). Von dieser Stelle ab wächst die dimensionslose Geschwindigkeit schwächer 
als linear mit 7. Sie entfernt sich mit steigendem ) immer weiter von der Geraden 7, um dann 
allmählich in den logarithmischen Verlauf überzugehen (etwa bei n = 100). 

Das Zurückbleiben der Geschwindigkeit gegenüber dem geradlinigen Verlauf n hat seinen 
Grund in dem Auftreten der Turbulenz. Den Einfluß der turbulenten Reibung erkennt man 
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T; 

e en all 
14 Zu 
das unter Verwendung von Gl. (12) in Bild 2 PR ist. Wie man aus beiden Abbildungen 
erkennt, ist eine nennenswerte turbulente Reibung erst ab 7 = 4 feststellbar, wo also das Zurück- 
bleiben der dimensionslosen Geschwindigkeit gegenüber ‘dem geradlinigen Verlauf n beginnt. 
Dort, wo der Übergang in die logarithmische Geschwindigkeitsverteilung erfolgt (also etwa 
bei n = 100), beträgt der Anteil der turbulenten Reibung an der gesamten Schubspannung 
über '98%. 

Die vorstehenden Betrachtungen haben selbstverständlich nur näherungsweise Gültigkeit, 
da wir die genaue Gesetzmäßigkeit des Abklingens der turbulenten Scheinreibung in Wandnähe 
nicht kennen. Der durch unsere 


% Formel (18) dargestellte Geschwin- 


digkeitsverlauf dürfte jedoch von 
der Wirklichkeit nicht wesentlich 
abweichen, da die Grenzbedingun- 
gen für die Austauschgröße sowie 
für die Geschwindigkeit sowohl an 
der Wand als auch in größerer Ent- 
fernung von der Wand richtig er- 
füllt wurden. 
Die experimentelle Prüfung 
des Wandgesetzes für kleine n ist 
kein einfaches Problem, da die Ge- 
oireichapdt schwindigkeitssonden (Stauröhr- 
chen, Hitzdrähte) in großer Wand- 
nähe und bei sehr kleinen Ge- 
schwindigkeiten Fehlanzeigen ge- 
ben, deren hinreichend genaue Er- 
mittlung große Schwierigkeiten 
bereitet. Solche bei kleinen n-Wer- 
ten durchgeführten Geschwindig- 
keitsmessungen, die der Verfasser 
15 | im Rahmen einer früheren Arbeit) 
a. 0 “ durchgeführt hat, sind in Bild 3 
Bild 3 . dargestellt. Wie man sieht, liegen 
die meisten Meßpunkte in der Nähe 
der gerechneten Kurve. Etwas größere Abweichungen treten lediglich bei den in sehr geringen 
Wandabständen durchgeführten Hitzdrahtmessungen auf und zwar in dem Sinne, daß die ge- 
messenen Geschwindigkeiten zu klein sind. 


ee 
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4. Das Mittengesetz 

Da im mittleren Teil eines Kanales oder Rohres alle Vorgänge symmetrisch zur Mittel- 
ebene bzw. Achse verlaufen, empfiehlt es sich, den Ursprung des Koordinatensystems in die 
Mittelebene bzw. in die Achse zu verlegen. Der Abstand von diesem Ursprung sei 2. z verläuft 
von r bis —r. Links von der Mitte im Abstande z = r möge sich die Wand befinden, für deren 
Nähe wir soeben die Geschwindigkeitsverteilung dargestellt haben. Es ist also 2= r—y bzw. 
yzr—12. 

Im mittleren Bereich kann die molekulare Zähigkeit gegenüber dem turbulenten Aus- 
tausch vernachlässigt werden. Führen wir noch —dz/r statt dy/r in unsere sung (5) 


ein, so erhalten wir 
u T 
\u To 


a) REN 
7 UN 
Über die dimensionslose Austauschgröße A/un, geben unsere Versuche Auskunft. Durch 


unmittelbare Messung des Differentialquotienten des Staudruckes mit Hilfe eines Doppelstau- 


5) H. Reichardt, Die Wärmeübertragung in turbulenten Reibungsschiehten. Z. angew. Math. Mal: 
20 (1940), S. 297. 


. (20). 
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£ rohres konnten wir die Austauschgröße vor allem für den mittleren Bereich ziemlich genau 
bestimmen (s. die Bilder 4 und 5)®). Diese Experimente lassen sich näherungsweise durch 
folgende Gleichung einheitlich darstellen: 


4 % z\?\/ 2\8 
un 3 08 (4) 2 (2) 
Er; + (054 : IK - SR UETN 
12 Führen wir diese Beziehung in (20) ein und be- 
achten, daß in ausgebildeter Strömung — EN = 
0 NEN 


ist, so folgt hieraus 


> ae d EA 
a 57 Rn —- konst (22). 
| +2) -6)) 


Die Integration ergibt 


Bin! a re TRIER ll: Ir N > 
un 5 [| tkonst . . . (28), 
& nz 
r 
070 r 
A_ 
Ar IH ——_ N _ 15 
V—0— Ui 
N N IL, 
005 nf? 
\ 
ou: 75, 2cm/s N 
® 19.8 cm/s \ 
8 46lcmjs are 
EN, \ 
r 
PER BERN SER) ER 0.25 gs 090 0.75 10 
Bild 4 Bild 5 


Die Konstante bestimmen wir für die Mittelebene bzw. Achse, wo - =(0 und v= u„ist: 
2 2 
FELL 
+25) 
ae MO) 
r 


Diese Formel bezieht sich auf den gesamten Strömungsquerschnitt, ausgenommen die 
unmittelbar an den Wänden liegenden dünnen Schichten, in denen die molekulare Reibung 
eine wesentliche Rolle spielt. Da in Gl. (24) kein Wandabstand, sondern nur der Abstand z 
‘von der Mitte eingeht, die die Symmetrieebene bzw. Symmetrieachse des Geschwindigkeits- 
‚profiles darstellt, haben wir es hier mit einem echten „Mittengesetz“ zu tun. 

Zur Prüfung der Formel (24) hat der Verfasser Geschwindigkeitsmessungen ausgeführt 
(s. Fußnote 6). Diese Messungen wurden in der durch Gl. (24) vorgeschriebenen Form 
aufgetragen. Wie Bild 6 zeigt (in dem statt des natürlichen Logarithmus der Zehnerlogarithmus 
erscheint) liegen die Meßpunkte auf bzw. in der Nähe einer Geraden mit der Neigung 5,75 (das 


. (4). 


entspricht einer Neigung von 4 — 2,5 bei Verwendung des natürlichen Logarithmus). An dem 
u — 0,9 
r 


unten eingezeichneten Maßstab für - ist zu erkennen, daß das Mittengesetz (24) über 
hinaus, also bis in die unmittelbare Nähe der Wände gilt. 


. . 6) In Bild 4 bedeutet P= ü/umar: — Die Messungen sind im mittleren Vertikalschnitt eines recht- 
eckigen Kanales durchgeführt worden, dessen Höhe 24,6 cm und dessen Breite 98 em beträgt. 


a - 
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Für kleine z/r läßt sich der Logarithmus in (24) durch eine Reihe ersetzen und man kann 


schreiben 
ee 
—1-- a ee 
PR Kun|ur\r ER 


Wir erhalten also für den mittleren Teil eine parabolische A des Geschwindigkeitsprofiles. 


In Bild’4, das die Geschwindigkeitsverteilung = für v= 15,2 cm/s zeigt, ist das - 


Max 


Parabelprofil der Formel (25) gestrichelt eingezeichnet. Wie man sieht, besteht zwischen 
dieser Parabel und der Kuppe des gemessenen Geschwindigkeitsprofils eine innige Berührung. 


1+2(2)° 
Eu e2 z a 
12 10 08 06 ” \r-(@) 02 v 


o u*12,2cm/s 
o _ 12cm/s 
®e  798cm/s 
o  220cm/s 
2 4ölcm/s 


Die parabolische Gestalt der Profilkuppe hat ihren Grund darin, daß sich der turbu- 
lente len! im Ha Teil des Strömungsquerschnittes kaum verändert 


(für x kleine — folgt aus (21) 


> = . Diesen Sachverhalt erkennt man auch unmittelbar ; € 
aus dem hraktisch linearen Verlauf ne Geschwindigkeitsgradienten im mittleren Strömungs- 


querschnitt. Denn da r— A n- ST, = ist, ist die Linearität von = ein Kennzeichen für 


die Konstanz der Austauschgröße. 


Es ist typisch für die turbulenten Profile, daß der NEE des Geschwindigkeitsgradienten 
mit wachsender Entfernung von der Mitte zunächst nachläßt, wie man aus Bild 4 ersieht, ; 


<0,7) liegen die Meßpunkte 


zeichneten linearen Verlaufes, den man in der Mitte beobachtet. Dieses Nachlassen der Profil- 
krümmung mit wachsendem z, das sich i in einem Flacherwerden des Profiles gegenüber der Para- _ 


bel äußert, hat seinen Fa in dem Anstieg der Austauschgröße (4 m ee r) (s. Bild 5). gr r IR 
Diese hat ihr Maximum bei un 0,5, wie bereits Nikuradse (s. oben) festgestellt. hat, Hier a 
2 > 


liegt daher auch das Minimum der Profilkrümmung 2 23 nn er} r N r e; 


Im mittleren Bereich (: unterhalb des gestrichelt einge- 


Wir wollen nun die errechnete Geschwindigkeitsverteilung (24) als Funlbion a) war. z “ 
abstandes darstellen und führen zu diesem Zweck den Abstand y von ey an > | 
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ein. Wegen y=r—z erhalten wir dann: 


(26). 


Wenn wir hier das zweite Glied der rechten Seite fortlassen, so erhalten wir das Mittengesetz 
in seiner bisherigen Form. 
‚Für die weitere Umformung von (26) definieren wir eine dimensionslose Zusatzgeschwindig- 


keit = durch folgende Gleichung: 


ee ER az 


* “ 
Führen wir ferner n, = — in (26) ein und bezeichnen den Wert von %; in der Mitte mit w;„, 
so folgt: 


U 1 %; Um 1 Um 
Fern RB EN 


Auf der rechten Seite dieser Gleichung stehen die Festwerte, die die einzelnen Glieder der linken 


- Seite in der Mitte des Rohres oder Kanales erreichen. Da die rechte Seite von (28) konstant ist, 


können wir auch schreiben 


7 1 u; 1 | 2 ( o5 z ) | 
BEN n Inn- Pe: - c= = In|n RN SICH (29), 
1+2[= 


| ) 


wobei c eine experimentell zu bestimmende Konstante bedeutet. 

Durch die Einführung des dimensionslosen Wandabstandes n in unsere Mittenformel (24) 
sind wir also zu einer Gleichung gelangt, die sich von dem bisherigen Wandgesetz um das zu- 
sätzliche Glied ,;/u* unterscheidet. Dieses Glied, das nur von der Ortskoordinate z/r bzw. y/r 
aber nicht von u* bzw. n abhängt, gibt die Abweichung der Geschwindigkeitsverteilung im 
Innern vom logarithmischen Profil wieder. Setzt man x = 0,4 in Gl. (27) ein, so erhält w,/u* 
den in Bild 5 dargestellten Verlauf. 

Wie Bild 5 zeigt, sind die Zusatzgeschwindigkeiten nur gering. Der bei y/r = 0,78 lie- 
gende Maximalwert von w;/u* beträgt 1,28. Von dort fällt «;/w* schnell ab bis auf den Wert 


1,01 in der Mitte. Dieser Verlauf gilt natürlich nur für den hier behandelten Fall der ausgebil- 


deten Strömung mit Druckabfall. Würde man eine entsprechende Funktion für die Grenzschicht 
an der ebenen Platte berechnen, so würde man, nach den Messungen von Schultz-Grunow 
zu urteilen?), auf etwa zweifach höhere Werte von u;/u* kommen. Diese höheren Zusatzgeschwin- 
digkeiten haben offenbar ihren Grund in dem Abklingen des turbulenten Austausches beim 
Übergang der Grenzschichtströmung in die äußere Potentialbewegung. 

Die eigenartige Form der Funktion w;,/w* erklärt die bekannten leicht gekrümmten Kurven- 
abschnitte, die die Meßpunkte im Inneren des Strömungsquerschnittes bei der halblogarith- 
mischen Auftragung beschreiben. Diese Abweichungen von der logarithmischen Geraden ver- 
u— U; \ 

ur 
hervor, in der die bereits in Abb. 6 wiedergegebenen Geschwindigkeitsmessungen in der geschil- 
derten Weise aufgetragen sind. Im Unterschied zu den “ -Meßpunkten sind die Punkte = 


als Funktion von logn darstellt. Das geht aus Bild 7 


4 > u 
schwinden, wenn man statt 
u 


ganz geradlinig ausgerichtet, woraus die Brauchbarkeit der Formel für w,/u* hervorgeht. 

Wie Bild 7 zeigt, liegen die Punkte — etwas unterhalb der gestrichelt eingezeich- 
neten Geraden, die die Neigung = — 5,75 hat und die Abszisse im Punkte ce = 5,5 schneidet. 
Auf eine genaue Ermittlung von c und x aus den vorliegenden Messungen müssen wir aber ver- 


?) F. Schultz-Grunow, Neues Reibungswiderstandsgesetz für glatte Platten. Lufo 17 (1940), 8. 239. 
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zichten, da die aus dem Druckgefälle bestimmten «*-Werte für diese Aufgabe zu ungenau sind). 
Wir können zunächst nur sagen, daß das Auftreten des Gliedes u;/u* grundsätzlich eine gewisse 
Voraussichtlich wird man den Wert x = 0,4 
beibehalten können und nur den c-Wert zu verringern brauchen. 

Auch in der neuen Auftragung ist, wie aus Bild 7 hervorgeht, noch eine gewisse Streuung 
der Meßpunkte vorhanden. Diese Schwankungen sind nur zum Teil durch Meßfehler bzw. durch 


Änderung der Konstanten c und x zur Folge hat. 


35 \25 


[e] 
® 
o 
a 


o u"72.2cm/s 


75,2 cm/s 
798 cm/s 
22.0. cm/s 
451 cm/s 


0 10 20 
i log 


Bild 7 


3.0 


Ungenauigkeiten bei der Mittel- 
wertsbildung bedingt. Wegen ge- 
wisser Mängel der Versuchsein- 
richtung (Störungen im Zustrom) 
waren die Geschwindigkeitsprofile 
zeitweise nicht völlig symmetrisch. 
Diese Störungen sind in die Mes- 


sungen mit eingegangen. Es ist 


aber noch etwas Anderes zu be- 
denken. 

Unsere in Bild 5 dargestell- 
ten Messungen von A/un, zeigen 
eine gewisse Abhängigkeit dieser 
Funktion von der Reynoldsschen 
Zahl bzw. von n,. Wenn diese Be- 


obachtung wirklich reell sein sollte, 


so wäre unsere einfache Glei- 
chung (21) (deren rechte Seite 
von n, unabhängig ist) für eine 
genaue Darstellung der turbulen- 
ten Profile grundsätzlich unzurei- 
chend. Diese prinzipiell außer- 
ordentlich wichtige Frage dürfte 
jedoch in diesem Zusammenhang 


keine Rolle spielen, da selbst durch eine so erhebliche Veränderung der Theorie, wie sie unser 
Ansatz (21) zur Folge hatte, nur die Tatsache zur Darstellung gelangen konnte, daß die Ab- 
weichungen von der logarithmischen Geschwindigkeitsverteilung im mittleren Teil des Strö- 


mungsquerschnittes nur gering sind. 


Wir haben nun noch die wichtige Frage zu beantworten, weshalb die logarithmische Form 
des Geschwindigkeitsprofiles so dominiert, obwohl das logarithmische Gesetz voraussetzungs- 


‚gemäß nur für die unmittelbare Nähe der Wand gelten sollte. Dieser merkwürdige Tatbestand 


hat seinen Grund darin, daß in der alten Theorie zwei Faktoren vernachlässigt worden sind, 
die, wenn man sie berücksichtigt, sich zum Teil gegenseitig aufheben. Der eine Faktor ist der 
Abfall der Schubspannung mit dem Wandabstand, der andere die Abweichung der Austausch- 
größe von dem angenommenen linearen Anstieg an der Wand. 

Um das zu an SR wir in unserer Grundgleichung (5) (in der 1 gegenüber A/u 


zu vernachlässigen ist) _ ——1—- (anstelle des früheren Wertes 1), sowie 
To 
u EL z . (80) 
AN r r 


anstelle des früher benutzten Ausdruckes AR, 


“( 
: 


Y 


Wir erhalten dann 


are 


Da hier Ke= 5 gekürzt werden kann, besagt diese Gleichung, daß die logarithmische Ge- 


schwindigkeitsverteilung bis zur Mitte des Strömungsquerschnittes exakt gelten würde, wenn für 
die Austauschgröße die Funktion (30) Gültigkeit hätte. Jedenfalls können die Abweichungen der 


5) Die Frage, ob und in welehem Ausmaß die Konstanten x und c von der Forn des Strömungsquer- r 
schnittes abhängen, kann nur durch weitere Versuche geklärt werden. 


tatsächlichen Geschwindigkeitsverteilung vom logarithmischen Profil auf die Abweichungen der 
gemessenen Werte A/u n, von der durch Formel (30) beschriebenen Funktion zurückgeführt werden. 


N nenn. Reichardt, Vollständige Darstellung der turbulenten Geschwindigkeitsverteilung 219 


Um diese Unterschiede beurteilen zu können, haben wir den Funktionsverlauf von (30) 
ebenfalls in Bild 5 eingetragen (siehe die gestrichelt gezeichnete Parabel). Wie man sieht, gibt 
diese Parabel die Tendenz des Verlaufes der Austauschgröße wesentlich besser wieder als die 
ständig steigende Gerade xy/r. Trotz teilweise erheblicher Abweichungen von den Meßpunkten 
stellt der parabolische Verlauf in Wandnähe sogar eine verhältnismäßig günstige Annäherung 
an die Wirklichkeit dar (etwa bis $ — 0,15reichend). In diesem Wandabstand hat die Strömung 
aber bereits eine erhebliche Geschwindigkeit erreicht (u = 0,7 Uyar). Die Geschwindigkeits- 
entwicklung vollzieht sich also größtenteils in einem Gebiet, in dem die logarithmische Verteilung 
ziemlich genau gilt.) 

Die restliche Geschwindigkeitsentwicklung (0,7< n< ) dehnt sich über ein weites 
Gebiet aus (0,15 < r = ı) . Im mittleren Strömungsquerschnitt ist also das Geschwindigkeits- 


profil relativ flach. Je kleiner aber der Geschwindigkeitsgradient ist, desto geringer ist auch der 
Einfluß des Wertes der Austauschgröße auf den weiteren Geschwindigkeitsverlauf. Wir ver- 
stehen daher, daß die im Inneren auftretenden größeren Abweichungen der Meßwerte Alun, 
von dem durch Gl. (30) dargestellten Verlauf eine so geringe Wirkung haben. 

Wie wir bereits gezeigt haben, erfolgen die Abweichungen von der logarithmischen Ver- 
teilung im Sinne zusätzlicher Geschwindigkeiten. Die tatsächlichen Geschwindigkeiten sind 
deswegen höher als die Geschwindigkeiten des logarithmischen Profiles, weil der wirkliche tur- 
bulente Reibungskoeffizient über dem größten Teil des Strömungsquerschnittes geringer ist 
als der dem logarithmischen Gesetz zugrundeliegende Reibungskoeffizient der Gl. (30). Wie 


man aus Bild 5 ersieht, liegen die gemessenen Werte Zn bis zu einem Wandabstand S 078 
T: 

unterhalb der durch GI. (30) dargestellten Parabel. Die Zusatzgeschwindigkeit ,/u* hat daher 

bei z — 0,78 ihr Maximum. Sie nimmt dann bei noch größeren Wandabständen bis zur Mitte 


wieder ab, da der wirkliche Austauschkoeffizient in diesem Gebiet größer ist als der. durch 
Gl. 30) dargestellte Koeffizient des logarithmischen Profils. 


' 5. Formel für den Gesamtquersehnitt 

Im 3. Kapitel haben wir eine Formel abgeleitet (Gl. (18)), die die Geschwindigkeitsver- 
teilung bis in die laminare Wandschicht hinein darstellt. Diese Geschwindigkeitsverteilung 
geht im Inneren in das logarithmische Profil über. 

Im 4. Kapitel haben wir gezeigt, daß das logarithmische Profil im mittleren Teil des Strö- 
mungsquerschnittes nicht genau gilt und haben die dort auftretenden zusätzlichen Geschwindig- 
keiten «;/u* berechnet (Gl. (27)). Diese Zusatzgeschwindigkeiten sind so definiert, daß (statt «/u*) 
U—U; 

ur 


dem logarithmischen Gesetz gehorcht. 


Wir können daher den Gültigkeitsbereich unserer Wandformel (18) auf den gesamten 
Strömungsquerschnitt ausdehnen, indem wir in Gl. (18) = durch en ersetzen. Wir erhal- 


ten dann, wenn wir die noch etwas unsicheren Zahlenwerte der Konstanten %, n, und b sowie 


- den von der Form des Strömungsquerschnittes beeinflußten Wert der Größe , = ce — — Inx 
wieder ausscheiden, & 


uw 1 | KDE RR 177 RR 8 >») 32 
A = hirn)+a(|l e zu EN A A 
Führen wir in (32) noch ü,/ü* aus (27) ein, so erhalten wir als endgültige Formel für die 
turbulente Geschwindigkeitsverteilung des Gesamtquerschnittes: 
2 
> ( ri ) | 


U 1 
ar zn (1 rel) +4 


N 99 
ER a) A a ER 
M (83) 

Wie man sich leicht überzeugen kann, geht diese Gleichung für kleine Wandabstände 
[£ -1) in das Wandgesetz (18) über. Andererseits folgt aus (33) für große 7 die Querschnitts- 


formel (29). 
Eingegangen am 28. August 1950. 


nl 
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Spannungsfunktionen achsensymmetrischer Deformation 


in Zylinderkoordinaten 
Von A. Timpe in Berlin 


Anschließend an frühere Arbeiten wird das Problem der torsionsfreien, achsensymmeirischen Defor- 
malion von Umdrehungskörpern mittels einer einzigen Bipotentialfunktion behandelt. Die Methode wird 
auf Zylinderkoordinaten übertragen und auf verschiedene früher untersuchte typische Plattenprobleme zwecks 
Gewinnung eines einheitlichen Ausgangspunktes für allgemeinere Fälle angewendet. 


The author has trealed, in an earlier paper, the problem of the torsion-free, axially-symmetrical defor- 
mation of bodies of revolution, by aid of a single bi-potential function. The method is extended to eylindric 
co-ordinates, and applied to several typical cases of plate deformation, examined in the earlier paper, to 
get a uniform startingpoint for the investigalion of more general cases. 


En continuant ses communications precedentes l’auteur traite le probleme de la deformation sans torsion 
et symetrigue par rapport 4 l’axe d’un corps de revolution & l’aide d’une seule fonction bi-potentielle. La 
methode est &iendue pour le cas de coordonne&es eylindriques et appliquee aux plusieurs probl&mes Iypiques, 
examines plus töt, de la deformation de plaques, dans l’intention de se procurer un point de depart pour de 
‚problemes plus gen£rales. 

IIponosnkası cBom mpeskHme paboTsI aBTop necaenyer upoÖnemy OcecuMmMeTpnyHoH 1ehop- 
Mann Te BpalmeHnst bes KPyYeHHA IIPM IIOMOINH TOIBKO ONHOH ÖMnoTenumHasBHoK PyHK- 
IHM. TOT METOA PACIPOCTPaHAeTcH Ha NHNHHNPMYeCKUC KOOPAHHATEI U IPHMEHAeTCH K yuce 
PAHBıme WHCeC/ICHOBAHHEIM XAPAKTePHBIM CAIYYAaAM TEOPHH IIIACTHHOK C MEJIBIO HONIYEHUA 
HUCXONHOH TOYKU MIA DONIee OÖMMX CIIYAeB. 


Das Problem der torsionsfreien achsensymmetrischen Deformation von Umdrehungskör- 
pern wurde vom Verf. in dieser Zeitschrift!) unter gleichzeitiger Verwendung zweier erzeugender 
Funktionen, einer Potentialfunktion 9 und der dazu assoziierten Strömungsfunktion y in 
Zylinderkoordinaten behandelt. Andererseits wurde von ihm später?) der Lösungsschlüssel einer 
einzigen der Bipotentialgleichung genügenden erzeugenden Funktion F in Polarkoordinaten 
aufgezeigt, wobei sich die Möglichkeit der Inversion achsensymmetrischer Verzerrungszustände 
ergab. Für den letzteren Fall wurden bald danach?) von ihm die charakteristischen ‚‚Brücken- 
lösungen“ des Bipotentialproblems diskutiert und insbesondere auf Kugelschalenprobleme an- 
gewendet. Die vorliegende Arbeit beschäftigt sich mit der Übertragung des Verfahrens einer 
einzigen erzeugenden Bipotentialfunktion F auf den Fall von Zylinderkoordinaten r, ©, z. Der 
Lösungsschlüssel, zu dem sie gelangt, wird auf verschiedene in der älteren Arbeit berührte 
typische Plattenprobleme angewendet (,‚Hauptlösungen‘‘) zwecks Gewinnung eines einheit- 
lichen Ausgangspunktes für allgemeinere Fälle. 


Grundgleichungen 
Sind im achsensymmetrischen Fall «,, 0, u, die Verschiebungskomponenten i in Richtung», 9,2, 
DEN u, au 
also, = —,‚ee=—,E,= "2 die Dehnungen, A die kubische Dilation, ferner 0,, Ce, 0,, T,, die 


or Falls Fear: 


Komponenten des Spannungszustandes, so haben wir die für das Volumelement geltenden 


Gleichgewichtsbedingungen 
00, , 0,—0e : OT, D 090 
r TE 4 —(: Lee ARTE 2 — Le eh lenge re . 
or 5. r r 02 i or r r F 02 er 


Sie werden durch die zwischen den Spannungs- und den Formänderungskomponenten bestehen- 
den Beziehungen 


=AA+2u,; 0e=AAt mE; = AAt u: her: 2) 
i Eo E 
Ed k FR LS a ee 3 9 
mit A aa u 3l+o) übergeführt in die Grundgleichungen: 
oA  1—2006 oA :1—2010(rö) 
rt Teder Ve N 
IE OU, 
Hierin ist 26 = Te — = - die einzige von null verschiedene KoiaionskunpZEE — Offenbar 


‘ 


PN 1—20 _ Re ’ 2 
besitzen A=g und a: Be die Eigenschaften einer achsensymmetrischen Potential- 


1 
funktion und der dazu assoziierten Strömungsfunktion: Durch Gleichsetzen der beiden aus den 
E 
Grundgleichungen zu gewinnenden Werte von 2 2: ergibt sich 
r 
2, u n% ve 
v7? > Ri ee = . . * * * . . EM 
1) Bd, 4 (1924), S. 361—376, ®) Z. angew. Math. Mech. Bd. 28 (1048), $. 161166, 


°) Z. angew. Math. Mech. Bd. 30 (1950), S. 50-61. 
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die bekannte Tatsache, daß die kubische Dehnung der Potentialgleichung genügt. Andererseits 


erhält man durch Gleichsetzen der beiden aus den Grundgleichungen zu entnehmenden Werte 
0° : Ä 
von I P_die Beziehung: 


or d2 ey Loy, My 


a 
FFE A EBE 1 BR EN TR FR NL RR I er a 
d.h. die für eine achsensymmetrische Strömung charakteristische Gleichung. 

Die Lösung der Grundgleichungen (3) wird nunmehr (in Anlehnung an den 1948 gegebenen 
Algorithmus in räumlichen Polarkoordinaten) auf zwei einer erzeugenden Funktion F ent- 
sprechende Anteile gestützt: 


Erster Anteil: 0; u = lü-o o)W?F+ 2 = 


F\ 
r)' 


Ihm entspricht: 


1 or u,) 2 — | ROVER ISO. 
7 r or Er E MEERE 02 r ordz2l’ 
EBK OU, OU, 1+o] a a 
SE PEwFT a a ar Wr 
s a Er LEW ®L u KR 
Zweiter Anteil = — , — gg el ae 
Ihm entspricht: 
Br a DR be 
Ba E rdrdz’ 
Di ER LeIO: 0° F 1-oor un 2) 
a da Tore arm Lör\ or 028 
Durch Kopplung der beiden Anteile entsteht: 
Br et oV®F 1+0o1 FF. 
d=e RS un C5) ET 
RR [er oV:F ee 5 er = 
285=—c (1—2 era Tas: N Im te ’ 
Für TERMIN sich: 
Sem a 
E 
2 F 
05 ge nn gs ne er 


Mit diesen Werten von A und 2& wird die»erste der beiden Grundgleichungen (3) identisch 
befriedigt. Andererseits ergibt sich durch Einsetzen in die zweite dieser Grundgleichungen die 


‚Forderung: 92 [a F rusE r3 "&]=0 


ya 2 a 1—o Ge) TuHOR: "Fr 


| ‚oder | en N NET nee (6), 


er 10ov®F &(V?F) 
or? VEEROT, EEE 


d.h, ‘die erzeugende Funktion muß der Bipotentialgleichung genügen®). Wir 
stellen die Grundformeln, auf die man sich bei gegebener erzeugender Funktion F zu stützen hat, 


VUN 


noch einmal zusammen: 


Nr | ji OF 
m ne — 2 Be 
iR or d2 ar re Ar F+ SR er rar’ 
OERLE Ir Ye au, , ro U, OU,, U, , OU, , 
Kr: a er e 
£ re j Ba Hg a E Eo E N 
re re “a ESEL ETee TEE: 0) Bf NILye 
te In Sir: 
3) Zu demselben Geshais ne auf umständlicherem Wege A.E.H.Love, Lehrbuch der Elastizität 


\ 7 » EL % 
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Wir überzeugen uns zunächst, wie die elementaren Deformationsfälle gleichförmiger 
Zug- bzw. Druckbeanspruchung eines Rotationszylinders aus der erzeugenden Funktion 


Bm or2 4 28. are a SE (8) 
gewonnen werden. (V?F=4cz2-+6z, also Br ü 
ER Ä ß Korn v8 
u, a u,=(1—2 or + ae 20) Moat6e)t Foren 
6, 420; 0- =— 20; = ne —(4c+0)(l—20) +40; 
Meet there (AH OF E 
Eo E Eco E 
ee a 7 EL 6)—2 ; 
Do) ae a 
Eo Eo E 
Br: 20.7 (1-0) (1—2o) Se a N ea 
br. er Eo E 
% Soll Radialspannung verschwinden, so muß a oO 
m sein, also | 
% LER Mn =, s ä 
F a ae N Br Met. 
% Soll Vertikalzugspannung verschwinden, so muß 0(4c+6)+ (1 —20)(4c+6)+4c=0 
EN sein, also 
A? LER LTE Em ENG = 
| ee ee 


Mit dem letzteren Ergebnis ist die erste Hauptlösung für das achsensymmetrische Problem 
der Kreisplatte im Rahmen der Bipotentialfunktionen aufgezeigt, mit der man allgemein die 
in die Ebene der Platte fallenden, am Rande resultierenden Zug- bzw. Druckkräfte erfaßt. 
Für das entsprechende Anliegen bezüglich der Platte, an deren Rande gleichmäßig ver- 
teilte Biegungsmomente wirken, benötigen wir erzeugende Funktionen vierten 


Grades. 
1. Harmonische Spannungsfunktion 4. Grades: 
e; Kein UretsM. u eenee: (9). 
Mit 7? F, = Oergibt sich gemäß (7): ne 2723 in er2en NR 
225 & = 225 &,= —Az2; 2.07% „sin; 
fr Der ae a kPa FE aM 
2. Biharmonische Spannungsfunktion 4. Grades: 
Ku= (tr 20); V:Fun=7r—142; ; 


u 2r2; u,—= (1 —20)7 (rt — 222) +87? — 22; = 22; = 2r; 


Dee. 4=—3(1—20)2;5 9,=(2—Bo)r; (9) 

7 ae er P . a 
0, n el 280-+2) 2; vr 280 +2)2; 8 RR ERS: 
,=— Ti (8 —70)4z; Bee; ie (8_70)2r 


Man kann die beiden Spannungsfunktionen Fr, Fir so kombinieren, daß 0, und EN ‚gleich. 
zeitig verschwinden. Damit ergibt sich die zweite nz ars für die a N 
ng Fir 2 “ \ 
w=(1—0o)-l14r2; 7er Us 
=0g=l4E:;. s,-0; 


ET 
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Auf die Längeneinheit des Kreisrandes der Platte von der Dicke 2 h erhält man als resultierendes 
TERE 28 ; . : 
Biegungsmoment 3 Eh?, Wie der Ausdruck für u, zeigt, wölbt sich die Mittelfläche z = 0 der 


Kreisplatte, die nur durch gleichmäßig über den Rand verteilte Biegungsmomente beansprucht 
wird, zu einer parabolisch gekrümmten Fläche. 


Die Kombination F** stellt die zweite Hauptlösung für die achsensymmetrische De- 
formation der Kreisplatte 2—= -- h, r= a dar, mit der man allgemein die Randmomente erfaßt. 
Hat man irgendeine den Bedingungen für die Oberfläche 2= + h genügende Lösung, so kann 
man die Befriedigung der Bedingungen bezüglich der Randmomente durch Überlagerung einer 
zweiten Hauptlösung erzielen. 

Als Material für die Weiterführung des Problems des achsensymmetrisch belasteten Kreis- 
zylinders, für den auf der Mantelfläche r—= «a die Normalspannungen und die Tangential- 
spannungen als Potenzreihen in 2 gegeben sind, vermerken wir die den Erzeugenden fünften 
Grades entsprechenden Lösungen. 


1. Harmonische Erzeugende 5. Grades: 


AR 
= 82 — 1022 4 15279) ; u rt ars = — 10242 22 
rR 
= or; 0 Ir; = 602-307; A=0; y,=120r2); 
E 1 ) E | 1 
> et RR}, FRE N; EN 
0, em ar S O9 ner 2 Jr . 
= a0 22% - 4 60 
3 en »ı-j@r2); u leg r2 
2. Biharmonische Erzeugende 5. Grades: 
Pu=4 824922); V?F—=9r2—6?°; 
16 
W=—2r —r2; u, = (17 — 180)r?2 + er 120)@ 
; 16 
5 —br— 2; 8 = —2r — 2; &,= (17 — 180)? — 3 3 120 au: 
'A4= 9(1— 20) (r? — 22) ; y,„=4(8—9o)rz, 
| E 21. 
> = Ir]; 0 =7 180) 22] ; 


= 57-90)" + (180 — 10)2] ee dern, 


1 ie 
Besonders wichtig für das Problem der achsensymmetrisch beanspruchten Kreisplatte 
- sind die geeignet normierten Erzeugenden sechsten Grades. 


1. Harmonische Spannungsfunktion 6. Grades: 


Hh=42—30#r+ ann  w=802r—60r2; = —40.-+ 12027 —15rt ; 


$ &, = 802 — 180 122 ; = Hr = — 1608 -+ 24072; 
a Be 9 | = 480 2 — 12078 ; 
RR RR E En 
Bu: TE =; — 180122) ; = 80; 


E 
h 0,= En ee 160 2? + 240 7?) ; T%, = I+e (240 2? r re 60r?) . 
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P2 


2. Biharmonische Spannungsfunktion 6. Grades: 


3 3 R, 4, . 
Fir Re s 


u,— — 12r?2 + 32729; 


V?Fı=3r— 247224828; 


u, = (1 — 20) (3r? — 24 r?2? + 824) + 241722? — 164; 


= —86r2 +22; = —12r2 4322; ,—=(1—20)(—481?2+ 322%) -- 48722 — 642° ; 


A=(1—20) (—48r°2+ 322); 


Eo E 


o (— 48 1?2 + 322°) - 


„= %(l+o)r2 — 2407; 


(— 36r?2 + 322°) ; 


"140 1-40 
Eo WER 
= —- „2 223) - er y2 9,23 . 
O9 age 4872-43 ET 12r?2 + 322°) ; 
16E 
le 
= we (1-+0)2?—30r?7] ; 


B 
u; = — [48 (1 + o)r2®? — 12072]. 


r 


Entscheidende Bedeutung für die durch Kombination der beiden Spannungsfunktionen 6. Grades 

zu gewinnenden Belastungsfälle hat nun die folgende Tatsache: Man kann die in den Aus- 
drücken für o, und r,, vorkommenden mit »?z bzw. r? proportionalen Glieder gleichzeitig zum 
Verschwinden bringen, wenn man die Spannungsfunktion # = Fr — 5 Fir zugrunde legt. Damit 
reduziert sich o, auf ein mit 2? und z,, auf ein mit r2? proportionales Glied. Sodann kann man Q 
diejenige aus den Spannungsfunktionen vierten Grades zu bildende Kombination heranziehen, 

für die 0, ,0, gleichzeitig verschwinden, während sich o, von der Form const.zund T,, von der Form 

const. r ergibt. Durch Überlagerung eines solchen Spannungszustandes über die Kombination . 
F=F}—5 Fır kann man olfenbar 7 für 2= + hzum Verschwinden bringen. Weiterhin kann man 
durch Überlagerung einer gleichmäßigen Zugspannung in der z-Richtung die Größe o, für 

z2—= — h annullieren. Wenn man schließlich noch das Biegungsmoment am Rande durch Hin- 
zufügen einer zweiten Hauptlösung beseitigt und durch Überlagerung eines konstanten Gliedes 


in u, das Verschwinden von u, für z=(0, r = 


a herbeiführt, so ergibt sich für die durch 


gleichförmige Last p auf z=h beanspruchte, in r=a frei aufliegende Kreis- 


platte folgender Deformationszustand: 
ur ateer 


Ur, 


Ban» 


1 


3 3 
E Ka.ht 3 EHE (++ a +9) Rt)! ; 


a a IE 18 Faro ,.| 


STAR 1— 0? I 


+ tes 
tert); en: | 
aut t+erte ee): 
eat ca E 
Ringegangen am 25. September 1950. 3 
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KLEINE MITTEILUNGEN. 


Zum Zirkeltest der Integrieranlagen. 

Meines Wissens existiert noch keine Theorie über die 
Fehler, die bei Integrieranlagen auftreten können. Um 
die Präzision von Reibradintegratoren zu untersuchen, 
werden diese sowohl einzeln wie auch in einfachen 
Schaltungen beobachtet. Insbesondere benutzt man 
bei vorwiegend qualitativen Untersuchungen eine 
Schaltung für die Differentialgleichung der harmoni- 


schen Schwingung, wie sie in Bild1 dargestellt ist. 


Es handelt sich um zwei Integratoren J, und J, mit 
Reibscheiben und Reibrädern. Es messe y; den Ab- 


2; den Drehwinkel des Reibrades beim Integrator J;. 


Y 


stand des Reibrades vom Reibscheibenmittelpunkt und 


Beide Reibscheiben seien mit demselben unabhängigen 

Antrieb verbunden; es sei @ der Drehwinkel der Reib- 
scheiben der beiden Integratoren. Bei passender Wahl 
der Maßeinheiten gilt . \ % 


da yıdaz = 1,2 u , j (51 ” R 


als Die ee der beiden Integratoren. 
‘ Nun ist die Schaltung in Abb. 1 so Desohatten, a 


r 
ot 


Na-ai- han... 
Bi ' Dr on 


& 
Pa Fr er 


Fi > > . A D 
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güt. Aus (1) und (2) folgt 


 rai=0; 


ty) 


Bild 1. 


Danach gilt allgemein 
2, = Asinz-+- Boos®; 


mit Konstanten A und 2, die sich aus den Abständen 
y, und y, der Reibräder für einen bestimmten Winkel x, 
etwa © = (0, berechnen. Werden z, und z, als karte- 
sische Koordinaten übereinander aufgetragen, so ent- 
steht ein Kreis, während die Schaltung angetrieben 
wird. Von einer guten Integrieranlage verlangt man, 
daß die gezeichnete Kurve sich schließt, wenn die Reib- 
scheiben um #—= 2 gedreht werden, Tatsächlich ist 
es möglich, Kreise von 40—50 om Durchmesser zu 
zeichnen, bei denen Anfang und Ende ineinander über- 
gehen und erst nach etwa 5—6 Umläufen Abweichungen 
vom ersten Kreise auffallen. 


Im folgendensolldargelegt werden 
daßals Folge von Getriebelose (back- 
lash) Spiralen geschrieben werden, 
die sioh von Umlauf zu Umlauf ver- 
breitern. 


$1. y, und y, sind nichtlineare Funktionen 
von 2, bzw. 2.. 

Wir untersuchen zunächst den Fall, daß Instrumen- 

tenfehler existieren, die die Relationen (2) verfälschen, 


jedoch nach wie vor y, als eindeutige Funktion von z,, 


y, als eindeutige Funktion von z, erscheinen lassen, Wir 


schreiben 
Yu -3— fl) = Fr (2’) 
sa At) vd) 
und setzen voraus, daß 9 (#) und y(#) für — o < x <® 
» stetig definiert sind, monoton zunehmen, und daß die 
Ungleichungen 
J@als<sM, lk@lsM ..0) 


mit einer gewissen Konstanten M bestehen. 
Zur Abkürzung setzen wir 


u. & 
[ v@)de = Fe); ward 09 (4). 


Dann folgt aus (1), (2°) und (4) 
i de, _ dz, 


Ten Fri 
dt, , dz, }- 
u’ Wtrz, et @)=0 
a und somit n 
Ri I EBEN H (2,25) >= F'(z,) -F G(z,) = (= oonst, (5). 
PD Wegen (3) gilt i 


I Ga)=-+ ©; 


Ie| + 


TEN 
a 2 % lim Fa)= +», 
N re ei 


Pi “ ferner besitzen F' (x) und. @ (2) N er der Monotonie von 
p und, Beben toner um. Nach links und 
nach reohts vom Minimum wachsen F(z) und G(x) 
nu »n gegen Unendlich, Auf die (2), 2,)-Ebene 

f\ 4 FL. N v 


ACH 
DE AND 
Pa M 


1.0 
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2= Bsinz — Aoose (4) 
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bezogen, existiert daher ein wohlbestimmtes Minimum i 
von H(z2,,2,), und längs jedem von der Minimalstelle 
ausgehenden Halbstrahl wächst H (2,, 2,) monoton gegen a 
Unendlich, Für C>Min / definiert daher die Gl. (5) ‘a 
eine geschlossene, doppelpunktfreie Kurve der (2, 2)- 
Ebene. Mit anderen Worten: | 
Der Zirkeltest ergibt geschlossene j 
Kurven. 

In (2’) gehen im wesentlichen die Fehler von Leit- 

spindeln ein, die zur Verschiebung der Reibräder dienen. 


$2., Der Einfluß der Getriebelose, 

Außer den zu (2’) führenden Fehlern, die eine Ab- 
weichung von den Relationen (2) verursachen, kommen 
im wesentlichen nur noch die Losefehler der Getriebe- 
leitungen in Frage. Der Losefehler läßt sich an einem 
einfachen mechanischen Modell leicht erklären. In 
Bild 2 bedeuten 1 und 2 zwei in Pfeilriehtung hin- und 
herbewegbare Teile, von denen das Teill efne Kerbe, 
das Teil2 eine in die Kerbe hineinragende Nase auf- 
weist. Wenn das Teill bewegt wird, nimmt es das 
Teil 2 mit; der Zusammenhang zwischen dem Weg «(t) 
des Teiles 1 und dem Weg v(t) des Teiles 2 (= Zeit) 
ist nicht linear, wenn die Nase die Kerbe nisht völlig 
ausfüllt. Es sei die Kerbe um die Länge 2a breiter als 
die Nase. Ferner sei 


ut)=Asintz; A>2u; (>20) (6), R 


und u = », wenn die Nase sich genau in der Mitte der 
Kerbe befindet, 
2 


Bild 2, 


Dann nimmt v () unter Vernachlässigung der Massen- 
trägheit des Teiles 2 den aus Bild 3 ersichtlichen Ver- 
lauf ©. Es ist 


v=u—«a zumindest für are sin = st<nj2, 


für 5 <t<h-+tnß, 


wobei 4 durch die sh <a; 


A sin (3 -F 1) = A — 2a definiert ist. Ferner ist 


Gleichungen 


v= 


u+ra fürn? Hh<ı<3n/2 =. 


v=-Ata rt <ı< rn 
v= A—d für S° + W"<siı<sön/2 


15 


A en 
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und so fort. Von gewissen Umgebungen der Extrem- 
werte von w(t) abgesehen, gilt also 


v(t) = ult) — asgnü (tl). (7) 
Es ist nützlich, das Integral 
ve han 
5 . 20 
= («(t) —_ v(t)) -sw(t)dt für © > aresin 1 
® 

. zu berechnen; diejenigen Intervalle, in denen v kon- 
stant ist, liefern keinen Beitrag zu J. Die anderen er- 


geben 
J=4a(A—a)!) (8) 


Dies vorausgeschickt, nehmen wir an, daß die Be- 
ziehung (2) durch Getriebelose verfälscht werden, Wir 
schreiben ‚jetzt 


"2 
y-5+42, = 42. . . (2), 

wobei 
Ariane, AR 0 SEh 2 aD) 


die Losefehler bis auf gewisse Umgebungen der Ex- 
tremwerte von z, und z, darstellen. 


Dabei kennzeichnet a,> 0 die Lose in der Getriebe- 
leitung vom Ausgang 2, zum Eingang Y,, a,>0 die 
. Lose der Getriebeleitung vom Ausgang 2, zum Eingang 
Y. Auch in den Phasen des Zirkeltests, in denen (9) 
nicht gilt, ist aber 


AN) 

erfüllt. 
Wenn x, der Drehwinkel der beiden Reibscheiben, 
monoton wächst, so kann man den Losefehler auf ® 


statt auf die Zeit i beziehen. Welchen Verlauf auch 
immer die Losefehler A2,(#) und Az,(z) nehmen mö- 


1) Diese Formelist noch weitgehender Verallgemeinerung fähig. 
Man bemerkt z. B., daß gewisse Variationen von u = A sin t für 
© <t<&-+2n das Integral nicht ändern, 
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gen, stets gilt 


008 (x -t)sin (x EN | 
Anl) Ana) ® 


17 
2, = Asin®+ B608% + 
6, (1) 
sin (@ —1)— 008 (8 -1) 


Az, (tl) Azs(t) a 


29 Bsin® -Ao00s#+ 


0 

mit geeigneten Konstanten A und B als allgemeine Lö- 
sung der Gleichungen (1) und (2’’), falls Az, und Az, als 
unabhängige Störfunktionen angesehen werden. Aus 
(10) und (11) folgt a 

2, —Asinae— Boosa| <ya3 + a2 © N (19) 

2, — Bsin x-+4A c00sa|< ya} +a}:® 
für #>0. Solange demnach 

ea Ha)<A+B=ÜC!; 0>0..(13) 

ist, können z, und z, in erster Näherung als Sinus- 
schwingungen angesehen werden. Alsdann läßt sich 
die Formel (8) anwenden, die auf eine Sinusschwingung 
Bezug nimmt. Es gilt also näherungsweise 
s+2nn 


H) Az; ();de zAna(C —a)zAna,l (14), 


B7 
wobei n eine ganze Zahl ist. Andererseits folgt aus (1) © 
und (2) 


d f : 
vn (4+2)= 22,42, 422542, . . (15). 
Zusammen liefern (14) und (15) die Formel 
vo r2rmn y 
(?2+2) | =8n(a, +0)0 (16). 
% 


Der Zirkeltest liefert daher, wie eingangs behauptet, 
eine sich verbreiternde Spirale, die in erster Näherung 
als Kreis mit dem Radius C angesehen werden kann, 
Bemerkenswert ist der Umstand, daß die Losefehler 
im gleichen Sinne auf eine Vergrößerung der Kreis- 
radien mit jedem Umlauf um den konstanten Betrag 
4 (a, -+ a,) hinwirken. Durch Beobachtung einer grö- 
Beren Zahl von Umläufen ist daher eine genaue Mes- 
sung von (a, + a,) möglich. 
Minden, Hans Bückner, 


BUCHBESPRECHUNGEN 


Federhofer, K. (o. Prof. a. d. Techn. Hochschule 
Graz, Prüfungs- und bungsauf- 
gaben aus der Mechanik des Punktes 
und des starren Körpers. In drei Teilen. 
‚: Teil: Statik. 165 Aufgaben nebst Lösungen. V -+- 
130 5. mit 243 Abb. Wien 1950. Springer-Verlag. 
Preis 9,60 DM. 


Das Studium der Mechanik erfordert außer der 
reinen Beschäftigung mit den theoretischen Grund- 
lagen vor allem die selbständige Bearbeitung von An- 
wendungsbeispielen. Eine hierfür geeignete Aufgaben- 
sammlung wurde im Jahre 1907 vom damaligen Grazer 
Mechanik-Professor F. Wittenbauer herausge- 
geben und erlebte im Jahre 1929 in der Bearbeitung 
von Th: Pöschlihre 6. Auflage, ist jedoch seit ge- 
raumer Zeit vergriffen. Die Aufgabensammlung von 
L.Fö bp l kam der regen Nachfrage bereits entgegen. 
Nunmehr legt der Nachfolger Wittenbauers in Graz, 
K. Federhofer eine eigene Aufgabensammlung 
vor, die auf dem Boden einer ca. 30-jährigen Lehr- 
tätigkeit erwachsen ist und sich in folgende Gebiete 
gliedert: Ebene Kraftsysteme und deren Gleichge- 
wicht, — Schwerpunkte ebener Flächen, — Ebene 
Fachwerke, — Biegemomente, Quer- und Längskräfte 
gerader Träger, — Dreigelenkbogen, — Raumkraft- 
system, — Seil- und Kettenlinien, — Stabilität des 
Gleichgewichts. Den weitaus größten Raum nehmen 
Beispiele der ebenen Statik ein, wobei auch die Er- 


‚der Spannungsverteilung in Konstruktionsteilen an 


mittlung der Verteilung sämtlicher Schnittgrößen 
(Biegemoment, Längs- und Querkraft) erörtert wird, 
wenn auch diese Größen erst in der Festigkeitslehre 
ihre Bedeutung erlangen. Bei Besprechung einzelner 
Aufgaben macht Verf. zuweilen auch von der Vektor- 
rechnung und anderen Hilfsmitteln (projektive Geo- 
metrie, Mayor-Mises-Abbildung) mit Vorteil Gebrauch. 
Bemerkenswert ist, daß der Abschnitt über Seillinien j 
auch Aufgaben über die Formbestimmung von Zylin- 
derschalen gleicher Festigkeit und über weit gespannte 
Kabel enthält. Das Verständnis wird durch ausrei- 
chenden Text und sorgfältig ausgeführte Abbildungen | 
erleichtert, so daß dem Studierenden mit dem vorlie- a 
genden Buche ein neues brauchbares Hilfsmittel zum 
Selbststudium in die Hand gegeben wird. 


Dresden. H.Neuber. 


Föppl, L. (o. Prof. a. d. Techn. Hochschule München) 
und E, Mönch (a.o. Prof. a. d. Univ. Tucuman), 
Praktische Spannungsoptik. VII + 
1628. mit 135 Abb. Berlin-Göttingen-Heidelberg 
1950. Springer-Verlag. Preis geb. 21,— DM. ; 

Die spannungsoptische Methode zur Untersuchung 


Modellen aus Kunstharz hat mehr und mehr an Be- 
deutung gewonnen. Besondere Verfahren der Modell- 
herstellung und -behandlung wurden entwickelt. Ds 
vorliegende Buch von L. Föpp1 und seinem Schüler 


a 


Fi 


ey 


ER 
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E. Mönch vermittelt die im Münchener Laborato- 
rium in den letzten Jahren gesammelten vielseitigen 
praktischen Erfahrungen und bildet so eine wichtige 
Ergänzung der bisherigen, in Deutschland erschiene- 
nen Literatur auf diesem Gebiete, welche im wesent- 
lichen nur die Grundlagen vermittelte. 


Dresden. H. Neuber. 


Dr.-Ing. Fritz Leonhardt und Dipl.-Ing. Wolthart 
Andrä, Die vereinfachte Trägerrost- 
berechnung. Mit 52 Hilfstafeln für die genaue 
und genäherte Berechnung von Rosten mit einem 
oder mehreren Querträgern, Angaben über Modell- 
untersuchungen und 7 Beispielen für die praktische 
Anwendung. 2. erweiterte Auflage. IV -+ 250 S. mit 
247 Abb. und Tab. Stuttgart 1950. Verlag Julius 
Hoffmann. Preis geb. 32,— DM. 

Die Verbindung zweier sich in der Regel recht- 
winklig kreuzender Scharen von Trägern hat in der 
jüngsten Zeit durch den Bau von Straßenbrücken 
und schwerbelasteten Trägerdecken besondere Be- 
deutung erhalten. Daher ist nach der Literatur wieder- 
holt versucht worden, die allgemein bekannten An- 
sätze zur statischen Berechnung der Schnittkräfte 
dieser hochgradig statisch unbestimmten Tragwerke 
für die üblichen Bauformen unter plausiblen An- 
nahmen zu zweckgebundenen Rechenvorschriften zu 
entwickeln und damit die Lösung der Aufgabe zu 
vereinfachen. Die Arbeit der beiden Verfasser ist eine 
Erweiterung der Anleitung für die vereinfachte Träger- 
rostberechnung von F. Leonhardt, dieim Jahre 1940 
erschienen ist. Die neue Auflage behandelt jedoch 
nicht nur den Trägerrost mit einem Querträger, son- 
dern auch mit mehreren Querträgern, wie dies in 
ähnlicher Weise von E. Metanund R. Schind- 
ler gezeigt worden ist. Der Trägerrost mit einem 
Querträger und beliebig einander zugeordneten und 
gestützten Hauptträgern kann ohne Schwierigkeiten 
nach den bekannten Ansätzen für die statische Unter- 
suchung des elastisch gestützten  durchlaufenden 
Trägers berechnet werden. Die Übertragung der 
Lösung auf den Rost mit mehreren Querträgern ist 
nur als Abschätzung zu bewerten, deren Brauchbar- 
keit für Vorberechnungen untersucht worden ist. Die 
genaue Lösung der Aufgabe beruht, kurz zusammen- 
gefaßt, auf der zweifachen Umordnung von Be- 
lastung und statisch unbestimmten Größen bei gleich- 
feldrigen Tragwerken, um die gegenseitige Abhängig- 
keit der Überzähligen zu beschränken. Die Grund- 
lagen sind ausführlich entwickelt, zu Rechenvorschrif- 
ten erweitert und durch zahlreiche Tateln und Ta- 
bellen zur unmittelbaren Berechnung von Träger- 
rosten erweitert worden. Den Abschluß der Arbeit 
bilden Angaben über die Untersuchung von Träger- 
rosten mit Modellen und Beispiele zur Anwendung 
der Theorie im Stahl- und Stahlbetonbau. Das Buch 
der beiden Verfasser ist eine wertvolle Bereicherung 

der Fachliteratur, das allen Fachleuten, die sich mit 
der allgemeinen Anordnung und mit der statischen 
Untersuchung von Trägerrosten befassen, gute Dienste 
leisten wird. 

Dresden. ‘K. Beyer. 

Dr.-Ing. habil. Ferdinand Schleicher (ehedem 

‘0. Prof. a. d. Techn. Hochschule Berlin), Taschen- 
buch für Bauingenieure mit Beiträgen zahl- 
reicher Fachleute. XXIII + 1942 S. mit 2403 Abb. 
Berichtigter Neudruck. Berlin-Göttingen-Heidelberg 
1949. Springer-Verlag. Preis geb. 36,— DM. 

Das Taschenbuch Für Bauingenieure ist der Fach- 
welt bereits durch die von M. Förster, weiland 


0. Professor a. d. Techn. Hochschule Dresden, heraus- 


gegebenen und in den Jahren 1911 bis 1928 erschiene- 


Buchbesprech ungen 
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schaftlich arbeitenden Bauingenieurs geworden. Die 
Erweiterung der Erkenntnis auf allen Teilgebieten 
der Bauingenieurwissenschaften hat zu einer voll- 
ständigen Neubearbeitung geführt, die F. Schlei - 
cher gemeinsam mit zahlreichen hervorragenden 
Fachleuten übertragen worden ist und zu einer 
wesentlichen Erweiterung und Vertiefung des Inhalts 
geführt hat. Das neue Taschenbuch ist im Jahre 1943 
erschienen, die Neuauflage vom Jahre 1949 ein be- 
richtigter Neudruck. Das Werk behandelt zunächst 
die Grundwissenschaften des Bauingenieurs, also 
Mathematik, technische Mechanik, Hydromechanik, 
Bodenmechanik, Festigkeitslehre und Elastizitäts- 
theorie. Das eigentliche Fachgebiet ist, um nur die 
wichtigsten Abschnitte zu nennen, durch die Bau- 
statik, die Baustofflehre, die Vermessungskunde, die 
Verkehrswirtschaft, den Straßenbau, das Eisenbahn- 
wesen mit Erd- und Tunnelbau, den Wasserbau mit 
der Wasserwirtschaft, den Grundbau, den See- und 
Hafenbau, den Städtebau und städtischen Tiefbau, 
den Massivbau und Stahlbetonbau, den Stahlbau 
und durch eine Einführung in die für den Bauingenieur 
besonders wichtigen Teilgebiete des Maschinenbaues 
vertreten. Alle Abschnitte tragen wissenschaftlichen 
Charakter und bieten durch ihre theoretische Fun- 
dierung und ihre Auswertung und Verbindung mit 
der Bauerfahrung weit mehr, als dies die bescheidene 
Bezeichnung ‚‚Taschenbuch‘ erwarten läßt. Es ist 
das Handbuch des Bauingenieurs, aus dem nicht nur 
der junge Student, sondern wohl alle, die sich diesem 
schönen Beruf gewidmet haben, bei allen Bauaufgaben 
mit Erfolg Rat und Hilfe holen werden. Dies wird 
durch die kurze, klare Zusammenfassung des Wesent- 
lichen, durch zahlreiche, gut ausgewählte Abbildungen, 
durch wertvolle Tabellen, Rechenvorschriften und 
Anwendungsbeispiele, vor allem aber durch die 
hervorragende Ausstattung unterstützt, die dem 
Taschenbuch für Bauingenieure auch in Zukunft die 
beherrschende Stellung in der Literatur der Bau- 
ingenieurwissenschaften sichern wird. 


Dresden. K. Beyer. 


Federhofer, K. (o. Prof. a. d. Techn. Hochschule 
Graz, Dynamik des Bogenträgers und 
Kreisringes. XII + 179 S. mit 35 Abb. und 
36 Zahlentafeln. Wien 1950. Springer-Verlag. Preis 
23,— DM. 

Nicht nur im Bauwesen, sondern auch im Maschinen- 
bau bilden die Bogenträger wichtige Bauelemente, die 
der wissenschaftlichen Festigkeitsforschung mannig- 
fache Anregungen gegeben haben. Die Statik der 
Bogenträger wurde in vielen Arbeiten des In- und Aus- 
landes analysiert und weitgehend für graphische Me- 
thoden zugänglich gemacht. Dagegen kann die Dyna- 
mik der Bogenträger noch nicht als abgeschlossenes 
Gebiet angesehen werden, so daß der Versuch des Verf. 
sehr zu begrüßen ist, die hier noch bestehende Lücke 
durch eine umfassende Darstellung auszufüllen. Dies 
gilt umsomehr, als er innerhalb der dynamischen Pro- 
bleme der Bogenträger und Kıreisringe eine weib- 
gehende Abrundung der wissenschaftlichen Durchdrin- 
gung erreichen konnte, wobei ihm vor allem eigene 
theoretische Arbeiten zustatten kamen, die zum Teil 
noch nicht an anderer Stelle veröffentlicht sind. 
Andererseits konnten die theoretischen Ergebnisse 
durch inzwischen bekannt gewordene Versuchs- 
berichte anderer Autoren bestätigt werden. Die Pra- 
xis wird das vorliegende Buch in vielfacher Hin- 
sicht begrüßen, da verwandte Probleme sowohl bei 
Schwingungen von Bogenbrücken wie auch von Lei- 
tungsrohren unter Innen- und Außendruck, ferner bei 
rotierenden elektrischen Maschinen und in ähnlichen 
Fällen auftreten. 


nen 5 Auflagen bekannt und durch seinen umfassenden Dresden. H. Neuber. 
Inhalt, zu einem wertvollen Hilfsmittel des wissen- 
Ile I Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen, 
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Nachrichten — Zuschriften an den Herausgeber 


Z. angew. Math, Mech, 
Bd.31 Nr. 7 Juli 1951 


NACHRICHTEN 


Tagung der Gesellschaft für angewandte 
Mathematik und Mechanik in Freiburg/Br. 


In der Zeit vom 28. bis 31. März 1950 fand die 
diesjährige Tagung der GaMM in Freiburg/Br. statt. 
Sie war von weit über 200 Teilnehmern besucht, unter 
anderem auch von zahlreichen Ausländern aus Eng- 
land, Frankreich, Holland, Italien, Österreich, Schwe- 
den, der Schweiz, den USA, Süd-Afrika usw. Nach 
den Ansprachen des Vorsitzenden der GaMM — Prof. 
Grammel, des Rektors der Universität Freiburg — 
Prof. Oelkers und des Oberbürgermeisters der Stadt 
Freiburg — Dr. Hoffmann, begann die eigentliche 
Tagung am Donnerstag Vormittag. Im ganzen wur- 
den etwa 75 Vorträge aus den Gebieten der prak- 
tischen Mathematik, der Elastomechanik und der 
Strömungsforschung gehalten. Die Themen dieser 
Vorträge hier anzugeben erübrigt sich, da eins der 
nächsten Hefte dieser Zeitschrift Vortragsauszüge 
und spätere Hefte einige Vorträge in extenso bringen 
werden. Der Donnerstag Vormittag und der Sonn- 
abend waren allgemein interessierenden Vorträgen 
gewidmet, während am Donnerstag Nachmittag und 
Freitag Parallel-Sitzungen abgehalten werden mußten. 
Am Freitag früh fand die Mitgliederversammlung der 
GaMM statt. In ihr wurde unter anderem mitgeteilt, 
daß die Mitgliederzahl im Laufe des letzten Jahres 
von 211 auf etwa 340 angewachsen sei, daß sie sich 
also um mehr als 50% erhöht habe. Der Mitglieder- 
beitrag wurde wieder auf DM 5,00 festgesetzt. 
Statt der aus den Vorstandsrat ausscheidenden 
Herren Pöschl und Weber wurden die Herren 
Hamel und Mettler gewäblt. Der Ort der 
nächsten Tagung konnte noch nicht bestimmt wer- 
den, doch soll die Tagung diesmal in Norddeutsch- 
land stattf inden,. Der durch, den Vorstand vollzogene 
Beitritt der GaMM zur JUTAM, die 1948 gegründet 
wurde, ein Band zwischen den einzelnen Fachorgani- 
sationen darstellt und deren Hauptzweck die Ver- 
anstaltung internationaler Tagungen ist, wurde ge- 
billigt, zumal die Aufnahme in diese Organisation 
einstimmig erfolgt war. Im September 1951 wird die 
JUTAM ein Symposium in der Riviera veranstalten, 
dessen Thema ‚Nichtlineare Schwingungen‘ ist. 


Für die Möglichkeit von Einzelaussprachen, für 
die ein besonderer Raum zur Verfügung stand und 
für Veranstaltungen zur Erholung war im reichlichen 
Maß gesorgt. Es fand am Mittwoch ein Empfangs- 
abend, am Donnerstageine Vorstellungim Großen Haus 
auf Einladung der Stadt Freiburg statt. Für Freitag 
Abend waren die Mitglieder von M. Pene (Commis- - 
saire pour le Land Bade) zum Cocktail im Deutsch- 
Französischen Club geladen. Am Sonnabend Nach- 
mittag fand ein Ausflug zum Schauinsland statt und 
am Sonntag machte dann noch ein Teil der Mitglieder 
eine Fahrt durch den Hochschwarzwald, ein anderer 
einen Ausflug nach Oberwolfach. Allgemein war man 
mit dem Verlauf der Tagung sehr zufrieden und es 
seiauch hier dem Organisator dieser Zusammenkunft — 
Herrn Prof. Dr. Görtler — für die mühevolle, sorg- 
fältige Vorbereitung dieser Tagung der herzliche Dank 
aller Teilnehmer ausgesprochen. Der bei dem letzten 
Abendessen insbesondere auch von den ausländischen 
Teilnehmern ausgesprochene Dank. wird ihm gezeigt 
haben, wie sehr man die große für die Vorbereitung 
der Tagung aufgewendete Mühe anerkannte. 


Willers. 


München: Am 26. April 1951 verstarb an den 
Folgen eines Verkehrsunfalles Geh. Rat Prof. Dr. 
Arnold Sommerfeld im 82. Lebensjahre. 


un 


Die Accademia dei Lincei hat im Dezember vorigen 
Jahres eine aus den Herren Oastelnuovo, 
Amaldi,Krall, Persico, Signorini und 
Tonolo bestehende Kommission beauftragt, die 
gesammelten Werke von Levi, Civita heraus- 
zugeben. Die Ausgabe wird etwa fünf bis sechs Bände 
umfassen. Der erste Band ist druckfertig und wird 
etwa 3000 Lire kosten. Um eine Entscheidung über 
die Höhe der Auflage treffen zu ‚können, bittet 
Herr Prof. A.Signorini- Rom (Universität) Inter- 
essenten sichanihn zu wenden. 


ZUSCHRIFTEN AN DEN HERAUSGEBER 


Druckfehlerberichtigung 


Zu Sauer: Dreidimensionale Probleme der 
Charakteristikentheorie partieller Differential- 
-gleichungen. Z. angew. Math. Mech. 30 (1950), 
8. 347— 356. 
Auf 8.354 findet sich die Bemerkung, daß die Diff.- 
Gl.(13) mit dem System der drei Diff.-Gleichungen: 


(8) Autı + .2.+ Ast G, 


Yu —- W= 


U u, =(0, 


äquivalent ist. Herr Pfeiffer hat mich freundlicher- 
weise darauf aufmerksam gemacht, daß diese Fest- 
stellung zu Mißverständnissen Anlaß geben kann. 
Natürlich sind die Diff.-Gl. (13) und das System (8) 
nicht in dem Sinne äquivalent, daß jede Lösung U, vd, W 

von ($) auch der Gleichung a, —v, —=(0 genügen und 
daher auf Lösungen der Diff.-Gl. (13) führen würde. 
Äquivalent sind jedoch die Lösungen von ($) und die 
Lösungen von (13) mit übereinstimmenden Anfangs- 
bedingungen, nämlich: Jede Lösung u, v, w von ($), 

deren . Anfangswerte u,, %y, W, auf einer raumartigen 
Ausgangsfläche (F) mit den Ableitungen 9, ®2, Ps 
einer Lösung 9 (x, y, 2) von (13) übereinstimmen, ist 
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in einer gewissen Umgebung von (F) identisch mit 
den Ableitungen dieser Lösung p und bestimmt daher 
diese bis auf eine additive Konstante. 


München. R. Sauer, 


Zu Buchholz: Besondere Reihenentwick- 
lungen für eine häufig vorkommende zwei- 
reihige Determinante mit Zylinderfunktionen 
und ihre Nullstellen. Z. angew. Math. Mech. 29 
(1949), S. 356— 367. 


S. 360: In Gl. (6a) ist im zweiten Faktor hinter dem 


Summenzeichen im Nenner z zu streichen. 


Zu Buchholz: Die Entwicklung willkür- 
licher Funktionen nach den Eigenfunktionen 
des Drehparabols. Z. angew. Math. Mech. 30 
(1950), S. 287—288. 

S. 288: In den GIn. (4) und (4a) muß ale Be 
tionsvariable n statt 7 gesetzt werden. 

In Gl. (5) lautet das Argument der M-HunEicuze 
unter dem Integral y statt @n. FR. 


In den Gin.(5) und (6) ist das Differential der, he} 


äußeren Integration dn durch d) zu ersetzen. 
Darmstadt. H. Buchholz. 


Isende an- 


en vor @ 


erkannt werden, anderenfalls wird das folgende Quartal noch geliefert. Verantwortlich für nn Anzeigenteil: Curt F. v ‚Schrei 


